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A view  of  instruments  in  the  cold  neutron  guide  hall  at  the 
NCNR  looking  back  towards  the  neutron  source.  The  blue 
biological  shielding  covering  cold  neutron  guide  NG-6  runs 
down  the  middle  of  the  picture.  In  the  foreground  are  two 
instruments  exemplifying  important  types  of  measurements 
on  materials  made  with  cold  neutrons:  on  static  atomic 
structures  (SANS,  red  tank  on  the  left),  and  on  atomic 
dynamics  (FCS,  below  the  yellow  structure  on  the  right). 
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DISCLAIMER 

Certain  commercial  equipment,  instruments,  or  materials  are  identified  in 
this  report  to  foster  understanding.  Such  identification  does  not  imply 
recommendation  or  endorsement  by  the  National  Institute  of  Standards 
and  Technology,  nor  does  it  imply  that  the  materials  or  equipment  identi- 
fied are  necessarily  the  best  available  for  the  purpose. 
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Foreword 


As  this  is  being  written,  we  have  just  begun  a three- 
month  shutdown  during  which  we  will,  among  other 
changes,  install  a next  generation  liquid  hydrogen  cold 
source,  which  will  increase  cold  neutron  intensities  by 
almost  a factor  of  two.  This  is  an  excellent  example  of 
the  necessary  operating  philosophy  for  any  facility  — 
improvement  and  change  are  the  only  constants.  During 
the  past  year,  we  have  had  many  changes  in  reactor 
systems,  in  instrumentation,  in  staff,  and  in  science,  all  of 
which  are  aimed  at  continuing  the  development  of  the 
NIST  Center  for  Neutron  Research  as  a first  rate  neutron 
facility. 

A new  cooling  tower,  which  is  specially  designed  to 
reduce  the  visible  vapor  plume  during  winter  months,  as 
well  as  to  provide  full  cooling  ability  for  the  next  twenty 
years,  has  been  constructed.  It  will  be  connected  into  the 
reactor  secondary  systems  during  the  shutdown  that  has 
just  started.  A new  analysis  of  major  reactor  systems  has 
been  initiated  using  the  most  modern  tools  as  part  of 
preparing  the  application  for  a license  extension  through 
2024.  Many  other  improvements,  ranging  from  instru- 
mentation renewal  to  better  provision  for  maintenance  to 
improved  electrical  distribution,  are  under  way  as  part  of 
the  preparation  for  re-licensing.  In  spite  of  the  shutdown 
that  just  started,  the  reactor  has  once  again  performed 
admirably,  operating  for  240  days  during  the  fiscal  year, 
meeting  our  scheduled  operating  time. 

The  three  new  inelastic  scattering  instruments,  which 
were  just  coming  on  line  last  year,  are  now  being  fully 
incorporated  into  the  user  program,  with  encouraging 
user  demand.  The  joint  National  Science  Foundation 
(NSF)/NIST  Center  for  High  Resolution  Neutron  Scatter- 
ing is  being  expanded  to  include  operation  of  these 
instruments  in  addition  to  the  prior  SANS  and  triple  axis 
instruments.  The  thermal  neutron  upgrade  program, 
which  includes  two  new  triple  axis  spectrometers,  an 
improved  11  Iter  analyzer  spectrometer,  and  other  instru- 
ments, is  advancing  well,  with  characterization  of 
monochromator  crystals  now  complete,  and  contracts 
awarded  for  the  manufacture  of  monochromator  drums 
and  focusing  devices.  A new  cold  neutron  triple  axis 


spectrometer  is  currently  entering  detailed  design  phase 
as  a joint  construction  project  with  Johns  Hopkins 
University  (partially  funded  by  the  NSF).  In  order  to 
create  more  research  space  in  the  guide  hall,  a new 
building  for  our  technical  support  activities  is  nearing 
completion.  In  response  to  many  recommendations,  we 
have  significantly  strengthened  our  sample  environment 
support  and  our  data  acquisition  and  analysis  computing 
effort,  adding  additional  staff  and  refocusing  efforts.  The 
number  of  research  participants  continues  to  grow, 
surpassing  all  of  our  expectations.  Overall,  we  are 
moving  forward  on  many  fronts  to  strengthen  our  re- 
search facility  operations  efforts,  in  order  to  better  serve 
our  customers. 

And.  as  always,  the  scientific  results  are  the  true 
measure  of  success,  and  the  highlights  that  form  the 
largest  portion  of  this  report  comprise  an  excellent  set  of 
examples  of  the  exciting  outputs  of  the  facility.  I hope 
that  you  enjoy  reading  them  as  much  as  I did. 
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Introduction  to  the  NIST  Center  for  Neutron 
Research  (NCNR) 


Modern  technological  society  is  dependent  upon 
increasingly  sophisticated  use  of  materials  with  attributes 
dictated  by  their  sub-microscopic  structural  and  dynami- 
cal properties.  A wide  range  of  techniques  provides 
knowledge  of  these  properties.  Scattering  probes  (for 
example:  x-rays,  light,  electrons,  neutrons)  are  arguably 
the  most  important.  Of  these  probes,  neutrons  are 
perhaps  least  familiar,  but  they  provide  important 
advantages  for  many  types  of  measurements. 

Modern  sources  provide  neutrons  moving  at  speeds 
comparable  to  those  of  atoms  at  room  temperature,  thus 
providing  the  ability  to  probe  dynamical  behavior. 
Neutrons  are  also  well  matched  to  probe  lengths  ranging 
from  the  distances  between  atoms  to  the  size  of  biologi- 
cal or  polymer  macromolecules.  Neutrons  are  sensitive  to 
the  magnetic  properties  of  atoms  and  molecules,  allowing 
study  of  the  underlying  magnetic  properties  of  materials. 
They  also  scatter  differently  from  normal  hydrogen  atoms 
than  they  do  from  heavy  hydrogen  (deuterium),  allowing 
selective  study  of  individual  regions  of  molecular 
systems.  Finally,  neutrons  interact  weakly  with  materials, 
providing  the  opportunity  to  study  samples  in  different 
environments  more  easily  (at  high  pressures,  in  shear,  in 
reaction  vessels,  etc.),  and  making  them  a non-destructive 
probe.  These  favorable  properties  are  offset  by  the 
weakness  of  the  best  neutron  sources  compared  to  x-ray 
or  electron  sources,  and  by  the  large  facilities  required  to 
produce  neutrons.  As  a result,  major  neutron  sources  are 
operated  as  national  user  facilities  to  which  researchers 
come  from  all  over  the  United  States  and  abroad  to 
perform  small-scale  science  using  the  special  measure- 
ment capabilities  provided. 

In  addition  to  scattering,  neutrons  can  be  captured  by 
nuclei  to  probe  the  atomic  composition  of  materials.  The 
subsequent  characteristic  radioactive  decays  provide 
“fingerprints”  for  many  atomic  nuclei,  allowing  studies 
of  environmental  samples  for  pollutants  (e.g.  heavy 
metals),  characterization  of  Standard  Reference  Materi- 
als, and  other  measurements.  There  are  important  areas  in 
physics  that  can  be  explored  by  measuring  fundamental 
neutron  behavior.  Examples  include  the  lifetime  of  the 
free  neutron,  important  for  the  theory  of  astrophysics;  the 


neutron  beta  decay  process,  that  provides  a stringent  test 
of  nuclear  theory;  and  the  effects  of  various  external 
influences  such  as  gravity  or  magnetic  fields  on  neutrons. 

The  NCNR’s  20  MW  NIST  Research  Reactor  pro- 
vides a national  user  community  with  facilities,  including 
the  nation's  only  internationally  competitive  cold  neutron 
facility,  for  all  of  the  above  types  of  measurements.  There 
are  about  35  stations  in  the  reactor  and  its  associated 
beams  that  can  provide  neutrons  for  experiments.  At  the 
present  time  28  of  these  are  in  active  use,  of  which  6 
provide  high  neutron  flux  positions  in  the  reactor  for 
irradiation,  and  22  are  beam  facilities.  The  following 
pages  show  a schematic  layout  of  the  facility.  More 
complete  descriptions  of  instruments  can  be  found  at 
http://www.ncnr.nist.gov. 

These  facilities  are  operated  both  to  serve  NIST 
mission  needs  and  as  a national  facility,  with  many 
different  modes  of  access.  Some  instrumentation  was 
built  years  ago,  and  is  not  suited  to  general  user  access; 
however,  time  is  available  for  collaborative  research. 
NIST  has  recently  built  new  instrumentation,  and  re- 
serves 1/3  of  available  time  for  mission  needs  with  the 
balance  available  to  general  users.  In  other  cases,  instru- 
mentation was  built  and  is  operated  by  Participating 
Research  Teams  (PRT);  PRT  members  have  access  to 
75  % of  available  time,  with  the  balance  available  to 
general  users.  Additionally,  NIST  and  the  National 
Science  Foundation  operate  the  Center  for  High  Resolu- 
tion Neutron  Scattering  at  the  NCNR,  with  one  thermal 
and  five  cold  neutron  instruments.  For  these  facilities, 
most  time  is  available  for  general  users.  While  most 
access  is  for  research,  with  results  that  are  freely  avail- 
able to  the  general  public,  proprietary  research  can  be 
performed  under  full  cost  recovery.  Each  year,  about 
1600  researchers  (persons  who  participated  in  experi- 
ments at  the  facility,  but  did  not  necessarily  come  here) 
from  all  areas  of  the  country,  from  industry,  academia, 
and  government  use  the  facility  for  measurements  not 
otherwise  possible.  The  research  covers  a broad  spectrum 
of  disciplines,  including  chemistry,  physics,  biology, 
materials  science,  and  engineering. 
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1.  A Cold  Neutron  Depth 
Profiling  instrument  (not 
shown)  for  quantitative 
profiling  of  subsurface 
impurities  currently  at  this 
site  will  be  moved  to  another 
position.  Shown  is  a 
proposed  Triple  Axis  Cold 
Neutron  Crystal  Spectrom- 
eter with  double  focusing 
monochromator  and  multiple 
crystal  analyzer/detectors 
that  can  be  flexibly  config- 
ured for  several  energies 
simultaneously  or  for  high 
throughput  at  one  energy. 

2.  BT-7  Triple  Axis  Spectrom- 

eter with  fixed  incident 

energy  for  measurements  of 
excitations  and  structure. 


3.  BT-8  Residual  Stress 
Diffractometer  optimized  for 
depth  profiling  of  residual 
stress  in  large  components 

4.  BT-9  Triple  Axis  Crystal 
Spectrometer  for  measure- 
ments of  excitations  and 
structure. 

5.  Thermal  Column  A very 
well-thermalized  beam  of 
neutrons  used  for  radiogra- 
phy, tomography,  dosimetry 
and  other  experiments 

6.  BT-1  Powder  Diffractometer 

with  32  detectors;  incident 
wavelengths  of  0 208  nm, 

0 154  nm,  and  0 159  nm, 
with  highest  resolution  of 
§d/d  = 8 x 10-1 2 * 4. 


7.  BT-2  Triple  Axis  Crystal 
Spectrometer  with 
polarized  beam  capability 
for  measurement  of 
magnetic  dynamics  and 
structure 

8.  BT-4  Filter  Analyzer 
Neutron  Spectrometer 

with  cooled  Be/Graphite 
filter  analyzer  for  chemical 
spectroscopy 

9.  BT-5  Perfect  Crystal 
Diffractometer  SANS 

small  angle  neutron 
scattering  instrument  for 
microstructure  on  the 
104  nm  length  scale 
sponsored  by  the  National 
Science  Foundation  and 
NIST,  part  of  the  Center 
for  High  Resolution 
Neutron  Scattering 
(CHRNS) 


10.  NG-7  Horizontal  Sample 
Reflectometer  allows 
reflectivity  measurements 
of  free  surfaces,  liquid 
vapor  interfaces,  as  well  as 
polymer  coatings 

11.  Neutron  Interferometry 
and  Optics  Station  with 
perfect  silicon  interferom- 
eter; vibration  isolation 
system  provides  excep- 
tional phase  stability  and 
fringe  visibility. 

12.  Spin  Polarized  Triple  Axis 
Spectrometer  (SPINS) 

using  cold  neutrons  with 
position  sensitive  detector 
capability  for  high  resolution 
studies  — part  of  CHRNS 
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13.  Spin  Echo  Spectrom- 
eter offering  neV  energy 
resolution,  based  upon 
Julich  design,  spon- 
sored by  NIST,  Julich 
and  ExxonMobil  — part 
of  CHRNS 

14.  Prompt  Gamma 
Activation  Analysis 

cold  neutron  fluxes  allow 
detection  limit  for  H of 
1 pg  to  10  pg  Focused 
beams  are  available  for 
profiling 

15.  NG-7  30  m SANS  for 

microstructure  measure- 
ments sponsored  by 
NIST,  ExxonMobil,  and 
the  University  of 
Minnesota, 


16.  Neutron  Physics  Station  a 

cold  neutron  beam  150  mm 
x 60  mm,  available  for 
fundamental  neutron  physics 
experiments. 

17.  Fermi  Chopper  TOF 
Spectrometer  a hybrid  time- 
of-flight  spectrometer  for 
inelastic  scattering  with 
incident  wavelengths 
between  0.23  nm  and 

0.61  nm  chosen  by  focusing 
pyrolytic  graphite  crystals  A 
simple  Fermi  chopper  pulses 
the  beam. 


18.  Disk  Chopper  TOF 
Spectrometer  versatile 
time-of-flight  spectrometer, 
with  beam  pulsing  and 
monochro-matization 
effected  by  7 disk 
choppers  Used  for  studies 
of  dynamics  in  condensed 
matter,  including  macromo- 
lecular  systems  — part  of 
CHRNS 

19.  NG-3  30  m SANS  for 

microstructure  measure- 
ments sponsored  by  the 
National  Science  Founda- 
tion and  NIST  — part  of 
CHRNS 


20.  Backscattering  Spectrom- 
eter: high  intensity  inelastic 
scattering  instrument  with 
energy  resolution  < 1 peV, 
for  studies  of  motion  in 
molecular  and  biological 
systems  — part  of  CHRNS 

21.  8 m SANS  for  polymer 
characterization,  sponsored 
by  NIST  Polymers  Division. 

22.  Vertical  Sample  Reflecto- 
meter  instrument  for 
measuring  reflectivities 
down  to  10'8  to  determine 
subsurface  structure,  with 
polarization  analysis 
capability 


NCNR  Images  2001 


Jeff  Nico  (center,  Ionizing  Radiation  Div.)  with  NRC  and  Harvard 
researchers  at  the  ultracold  neutron  trapping  experiment 
to  measure  the  neutron  lifetime. 


Andreas  Michels  (U.  Saarlandes)  and  Raghavan  Viswanath 
(Fz.  Karlsruhe)  at  the  NG-3  SANS  display. 


NCNR’s  John  Copley  describes  the  NIST  reactor  to  a visitor. 


Syracuse  University  Professor  Bruce  Hudson  (right)  with  his  team 
at  the  Filter  Analyzer  Neutron  Spectrometer. 


NCNR’s  Steve  Kline  injects  a sample  into  the  SANS  high 
pressure  cell. 
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Gathering  of  Eagles  -- 
in  search  of  the  elusive 
CDW. 


NCNR’s  Susan  Krueger  and 
Deborah  Kuzmanovic  (Bio- 
technology Div.)  collect  SANS 
data  on  molecules  that  can  be 
used  to  test  for  diseases  in 
humans. 


Jamie  Schulz  (University 
of  Maryland)  mounts  a 
sample  on  the  NG-1  re- 
flectometer. 


NCNR  Support  Staff  enjoy  a day  away  in 
Washington  to  visit  the  Dept,  of  Commerce 
and  the  White  House. 
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NCNR’s  Paul  Liposky,  Dave  Anderson  (FDA),  and  Liz  Mackey  (Neu- 
tron Methods  Group)  with  the  new  thermal  neutron  Prompt 
Gamma-ray  Activation  Analysis  Instrument  (PGAA). 


NCNR’s  Nick  Rosov  instructing  at  the  7th  annual  NCNR  Sum- 
mer workshop. 


Measurement  Of  The  Free  Neutron  Lifetime 


Free  neutrons  arc  unstable  and  decay  with  a lifetime 
of  about  15  minutes  into  a proton,  an  electron,  and  an 
antineutrino.  Accurate  determination  of  tins  natural 
lifetime  is  important  for  tests  of  the  Standard  Model  of 
Fundamental  Particles  and  Interactions,  as  well  for 
understanding  certain  aspects  of  cosmology  and  astro- 
physics. The  neutron  lifetime  influences  the  predictions 
of  the  theory  of  Big  Bang  nucleosynthesis  for  the  primor- 
dial helium  abundance  in  the  universe  and  the  number  of 
different  types  of  light  neutrinos.  Since  the  dominant 
uncertainty'  for  the  primordial  helium  abundance  comes 
from  the  accuracy  of  the  neutron  decay  rate  | I |.  im- 
proved neutron  decay  rate  measurements  are  needed  for 
sharpening  this  prediction. 

The  measurement  of  the  lifetime  of  free  neutrons  in  a 
cold  neutron  beam  was  one  of  the  first  experiments 
operating  in  the  NCNR  Guide  Hall  The  experiment 
occupied  the  NG-6  Physics  Station  alternatively  with 
other  fundamental  pin  sics  experiments  through  Decem- 
ber 2000.  A statistical  precision  of  one  part  per  thousand 
was  achieved  in  the  final  result.  An  overall  standard 
uncertainty  of  about  four  parts  per  thousand  has  been 
obtained,  and  that  uncertainty  may  be  eventually  reduced 
by  as  much  as  a factor  of  two  by  ongoing  developments 
in  characterization  of  the  neutron  counting  for  this 
experiment. 

Figure  1 shows  a schematic  diagram  of  the  apparatus. 
An  axial  5 Tesla  magnetic  field  provides  radial  contain- 
ment of  the  decay  protons,  while  axial  containment  is 
achieved  by  the  positive  electrostatic  charge  on  two  end 
electrodes,  called  the  mirror  and  the  gate.  This  combina- 
tion of  magnetic  and  electrostatic  fields  is  called  a 
"Penning  Trap.  The  neutron  detector  consists  of  a well- 
characterized  isotopic  target  of  "Li.  viewed  by  a set  of 
four  charged  particle  detectors  The  charged  particles 
from  neutron  reactions  in  the  isotopic  targets  arc  counted 
with  an  accuracy  of  better  than  one  part  per  thousand 
The  reaction  cross  sections  for  these  standard  cross 
sections  arc  known  to  about  0.25  %. 

When  a neutron  in  the  beam  decays  inside  the  trap- 
ping region,  the  recoil  proton  is  trapped.  Periodically,  the 
gate  electrode  is  lowered  to  ground  potential  to  allow  the 
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Fig.  1.  Schematic  drawing  of  the  NIST  Penning  trap  neutron  life- 
time experiment 

accumulated  trapped  proton(s)  to  exit  the  trap  and  be 
counted  by  the  proton  detector.  A bend  in  the  magnetic 
field  makes  it  possible  to  locate  the  proton  detector 
outside  of  the  incoming  neutron  beam.  The  trapped 
protons  have  a maximum  energy  of  less  than  800  eV  and 
arc  undetectable  at  those  low  energies,  but  by  holding  the 
proton  detector  at  a negative  high  voltage  of  the  order  of 
-30  kV.  the  protons  may  be  detected  with  a signal  well 
above  noise  lex  els.  The  signal-to-noisc  ratio  is  further 
enhanced  by  counting  for  only  a short  time  (less  than 
100  ms)  after  the  trap  is  opened  in  each  counting/ 
trapping  cycle.  Typical  trapping  times  are  of  the  order  of 
1 0 ms.  and  typical  decay  event  rates  are  of  the  order  of  a 
few  protons  per  second. 

Determination  of  the  lifetime  requires  either  an 
accurately  known  trap  length  or  a variable  trap  length 
with  accurately  known  differences  in  length  The  latter  is 
more  easily  realized  physically  and  is  achieved  in  this 
apparatus  by  the  precision  machining  of  the  segmented 
trap  structure.  The  hard-to-define  end  effects  do  not 
change  significantly  as  the  trap  length  is  varied  by  applying 
the  mirror  voltage  to  different  segments.  Figure  2 shows  a 
plot  of  the  measured  proton  rate  (after  normalization  to 
the  neutron  rate)  as  a function  of  trap  length;  the  slope  of 
this  line  is  inversely  proportional  to  the  neutron  lifetime. 

Since  the  time  when  some  of  the  present  staff  of  this 
experiment  participated  in  a closely  related  experiment 
1 2 1 at  the  Institut  Laue-Langevin.  a large  number  of 
changes  have  been  made  to  improve  the  accuracy. 

Perhaps  the  most  important  factor  w as  simply  having 
much  more  beam  time  for  evaluation  of  significant 
systematic  effects  and  accumulation  of  better  statistical 
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Fig.  2.  A plot  of  the  proton  count  rate  (normalized  to  the  neutron 
monitor  count  rate)  vs.  trap  length.  The  slope  of  the  line  is  in- 
versely proportional  to  the  measured  neutron  lifetime. 

precision.  Other  factors  include  Monte  Carlo-computed 
corrections  due  to  magnetic  field  inhomogeneities, 
improved  beam  profiles  and  proton  alignment,  better 
trap  voltage  stability  and  monitoring,  improved  defini- 
tion of  the  areal  mass  density  of  the  isotopic  targets,  use 
of  a new  trap  designed  to  minimize  instabilities,  and 
improved  analysis  methods. 


Backscattering  Fraction 

Fig.  3.  A linear  fit  of  the  measured  neutron  lifetime  at  varying 
values  of  the  detector  backscattering  fraction.  The  extrapola- 
tion to  zero  backscattering  gives  the  free  neutron  lifetime. 


Figure  3 shows  the  extrapolation  of  the  measured 
lifetime  as  a function  of  the  backscattering  fraction  for 
protons  from  the  detector  surface  to  zero  backscattering. 
The  preliminary  lifetime  value  obtained  in  this  experi- 
ment [3 1 is 

885.3  s±  4.0  s. 

Continuing  efforts  to  measure  the  neutron  count  rate 
are  underway  by  both  calorimetric  and  coincidence 
techniques,  which  should  reduce  the  final  uncertainty  by 
half.  They  will  not  only  improve  the  results  of  this 
experiment  but  also  provide  a more  accurate  and  direct 
calibration  of  the  NIST  Standard  Neutron  Source. 

Also  underway  at  NG-6  is  a new  experiment  to 
measure  the  neutron  lifetime  by  an  independent  method 
involving  magnetic  trapping  of  ultracold  neutrons.  This 
experiment  is  being  carried  out  in  collaboration  with  the 
Physics  Department  of  Harvard  University. 
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ETS-4:  Maxwell’s  Zeolite 


Scottish  physicist  James  Clerk  Maxwell  ( 1831-1879) 
discussed  a mythical  device  that  could  separate  mol- 
ecules by  temperature,  a concept  now  known  as 
Maxwell’s  demon.  Maxwell  used  this  idea  to  explore 
how  sorting  molecules  by  size  would  violate  the  second 
law  of  thermodynamics.  If  one  did  have  a Maxwell's 
demon,  a more  useful  task  for  it  would  likely  be  to 
separate  molecules  by  size.  This  would  allow  mixtures  to 
be  separated  into  pure  substances  and  would  facilitate 
removal  of  dangerous  contaminants  from  valuable 
materials.  In  some  cases,  molecular  separations  can  be 
performed  using  zeolites  and  similar  materials.  Zeolites 
are  composed  of  linked  tetrahedral  S i04  and  A104 
species,  arranged  in  rings  to  create  a framework  of 
molecule-sized  cavities.  Related  zeolitic  materials  may 
use  other  tetrahedral  atoms  to  build  the  framework. 
Molecules  can  access  these  cavities  or  pores,  under  two 
conditions:  if  they  can  pass  through  the  appropriate  rings 
in  the  framework  and  these  rings  are  not  blocked  by  other 
species  in  the  material,  such  as  cations.  If  the  ring  size  in 
a material  matches  the  needs  for  a particular  chemical 
separation,  the  material  may  be  appropriate.  Chemical 
modifications,  such  as  ion  exchange  or  changing  the 
framework  composition  can  sometimes  influence  the 

Room  Temp  150°C  200°C  250°C  300°C 

Fig.  1.  Size  of  the  eight-ring  pore  opening  in  ETS-4  as  a function 
of  dehydration  temperature. 

pore  size.  So,  occasionally  it  is  possible  to  modify  a 
zeolite  to  a particular  separation.  However,  this  does  not 
allow  what  Maxwell  might  envision,  where  the  pore  size 
can  be  adjusted  over  a wide  range  to  suit  the  desired 
process,  in  effect  a Maxwell's  demon  for  size  separa- 
tions. This  exact  capability  has  now  been  demonstrated  in 
material  ETS-4  [I  ].  ETS-4  can  be  processed  in  such  a 
way  so  that  the  pore  sizes  can  be  tuned  for  particular 
separations. 


ETS-4  is  a titanosi  1 icate  material  invented  by 
Engelhard  in  the  I980's  [2],  ETS-4  differs  from  most 
zeolitic  materials  in  a number  of  ways.  Chemically,  it  is 
an  oxide  of  Ti  and  Si,  while  the  majority  of  zeolitic 
materials  are  oxides  of  A!  and  Si.  Structurally,  the 
material  is  composed  of  both  tetrahedral  Ti04  and 
octahedral  TiO(  units,  as  well  as  SiOt  units.  Finally,  while 
most  zeolitic  materials  hydrate  and  dehydrate  reversibly 
or  are  destroyed  during  dehydration.  ETS-4.  at  least  in 
some  ion-exchanged  forms,  can  be  induced  to  dehydrate 
irreversibly  and  without  destruction  of  the  pores.  Dehy- 
dration reduces  the  average  pore  opening  size,  and  thus  it 
is  possible  to  tune  the  zeolite  to  accept  molecules  of  a 
particular  size  by  changing  level  of  dehydration.  The 
amount  of  dehydration  in  turn  is  dictated  by  the  tempera- 
ture to  which  the  material  is  exposed.  In  1999.  the 
Advanced  Technology  Program  funded  Engelhard  to 
exploit  this  effect,  which  Engelhard  calls  the  Molecular 
Gate,®  to  develop  methods  for  separating  oxygen  from 
air,  a process  with  useful  applications  for  prevention  of 
pollution  as  well  as  medical  and  other  industrial  applica- 
tions. 

The  mechanism  behind  the  Molecular  Gate®  effect 
was  demonstrated  crystallographically  using  data  from 
the  BT-1  neutron  powder  diffractometer.  Pores  in  ETS-4 
must  be  accessed  via  a ring  of  eight  Si  and  Ti  atoms,  each 
linked  by  oxygen  atoms.  This  ring  is  commonly  called  an 
eight-ring,  despite  the  fact  that  it  is  actually  composed  of 


van  der 

Waals 

opening 

(A) 

Dehydration  Temperature,  °C 

RT 

150 

200 

250 

300 

Di 

4.27 

3.97 

3.95 

3.94 

3.90 

d2 

4.43 

4.02 

4.09 

4.57 

4.57 

d3 

3.61 

3.28 

3.29 

3.27 

2.77 

Table  1.  Eight-ring  opening  as  a function  dehydration  tempera- 
ture. Distances  are  the  separation  between  van  der  Waals  radii 
between  opposite  05  (Dt),  01  (D2)  and  02  (D3)  atom  pairs  (see 
Fig.  2). 
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sixteen  atoms.  In  Fig.  1,  the  size  of  this  eight-ring  is 
shown  as  a function  of  dehydration  temperature.  The 
actual  size  of  molecules  that  can  be  admitted  via  this 
opening  are  dictated  by  the  van  der  Waals  radii  of  the  O 
atoms  in  the  eight-ring  and  the  van  dcr  Waals  dimensions 
of  the  molecule  to  be  admitted.  Since  accurate  determina- 
tion of  oxygen  atom  siting  is  needed,  and  since  ETS-4 
does  not  form  single  crystals,  neutron  powder  diffraction 


Sil  05  Sil 


Fig.  2.  Atom  and  pore-opening  distance  labeling  scheme  in  the 
eight-ring  opening  of  ETS-4. 

measurements  are  the  only  way  to  obtain  accurate 
measurements  of  the  ring  openings  as  a function  of 
dehydration  temperature.  Measurements  of  the  eight-ring 
opening  from  Rietveld  refinements  are  shown  in  Table  1 . 
relative  to  the  distances  diagrammed  in  Fig.  2. 

The  proof  that  dehydration  indeed  docs  change  the 
pore  access  is  demonstrated  by  adsorption  isotherms, 
shown  in  Fig.  3.  The  separation  shown  in  the  top  of  this 
figure,  between  O,  and  N show  s how  precisely  the  pore 
sizes  can  be  adjusted.  The  difference  in  van  der  Waals 
radius  between  oxygen  and  nitrogen  is  only  0.1  A. 


Fig.  3.  Adsorption  of  selected  molecules  as  a function  of  dehy- 
dration temperature. 
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Stresses  in  Electron  Beam  Weld  Joints  of 

Superalloys 


Single-crystal  turbine  blades  are  the  pinnacle  of  a 
decade-long  development  and  refinement  of  nickel-based 
superalloys.  By  eliminating  creep  along  grain  boundaries 
as  the  main  failure  mechanism  of  polycrystals  at  high 
temperatures,  the  use  of  single  crystal  turbine  blades  has 
allowed  a rise  in  temperature  in  the  "hot  zone’  by  more 
than  100  K.  thus  substantially  increasing  the  thermody- 
namic efficiency.  The  complete  turbine  blade,  including 
the  mounting  base  and  the  cooling  channels,  is  grown 
from  the  melt  in  a single  process  and  represents  a 
considerable  investment.  In  order  to  protect  this  invest- 
ment. electron  beam  welding  (see  Fig.  1 ) is  used  to  repair 
minor  structural  damage  of  failed  components. 

[022] . x 


microstructure 

P . 


s.#-. 


Fig.  1.  Test  specimen  of  a weld  joining  a Rene  N5  single  crystal 
and  a polycrystal.  The  three  different  regions,  polycrystal,  weld, 
and  single  crystal,  appear  in  different  colors  after  etching  which 
indicates  a different  chemical  composition.  The  microstructure 
shows  typical  features:  60  % cuboidal  / -precipitates  and 
unimodal  size  distribution. 

Usually  a polycrystalline  alloy  of  comparable  compo- 
sition is  chosen  as  a replacement  material.  The  typical 
grain  size  of  the  polycrystal  is  approximately  one  milli- 
meter which  is  about  the  size  of  the  incident  neutron 
beam.  Thus,  stress  measurements  arc  not  feasible  for  this 
side  of  the  weld. 

The  weld  process  creates  very  high  residual  stresses 
that  make  the  extended  weld  region  prone  to  cracking. 
Cracks  arc  sometimes  found  parallel  to  the  weld  espe- 


cially in  the  z-dircction.  so  that  an  additional  heat 
treatment  must  be  applied  for  stress  relief.  The  standard 
heat  treatment  is  2 h at  1423  K for  homogenizing  (solu- 
tion) and  a subsequent  aging  treatment.  In  this  investiga- 
tion we  compare  the  stresses  before  and  after  heat 
treatment  in  the  single  crystal  side  of  the  weld  joint. 

Compared  to  a measurement  on  a fine-grained 
polycrystal,  the  measurement  procedure  is  rather  difficult 
because  all  small  angle  grain  boundaries  (Fig.  2)  have  to 
be  brought  into  reflection  and  averaged  over  the  rocking 
angle  co  in  the  data  reduction. 


Fig.  2.  Q^-intensity  map  for  the  as-welded  sample  (reciprocal 
space)  plotted  vs.  26 and  the  rocking  angle  co.  Small  angle  grain 
boundaries  as  indicated  by  the  presence  of  several  maxima  are 
characteristic  for  superalloy  single  crystals.  The  existence  of  two 
maxima  in  20  (equivalent  to  two  d-spacings)  indicates  the  /-ma- 
trix and  the  /'-precipitates  whose  lattice  parameters  are  close  to 
each  other. 

Additionally,  both  the  welding  and  the  heat  treatment 
produce  position-dependent  chemical  gradients  that 
affect  the  unstressed  lattice  parameters.  Their  measure- 
ment would  have  required  cutting  the  specimen  into  little 
pieces  for  complete  stress  relief.  We  chose  an  alternative 
approach  of  applying  mechanical  boundary  conditions  at 
the  surface.  The  results  obtained  from  this  procedure  are 
shown  in  Fig.  3. 

The  effect  of  the  heat  treatment  is  dramatic.  Both  the 
strong  tensile  stresses  in  the  heat-affected  zone  (5  mm 
from  the  weld  center)  and  the  strong  gradients  at  the 
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interface  between  weld  and  heat  affected  zones  are 
removed.  Stress  relief  in  the  z-direction  is  especially 
important  because  it  removes  the  driving  force  for  crack 
opening  parallel  to  the  weld  in  the  heat-affected  zone. 


stresses:  z-direction 


distance  from  weld  middle  [mm] 


Fig.  3.  Residual  stresses  (th  before  and  after  the  heat  treatment 
for  two  different  depths  from  the  upper  surface.  AW  stands  for 
‘As  Welded’;  AS  stands  for  ‘Heat  Treated’.  The  uncertainties  are 
of  the  order  of  the  size  of  the  data  points. 

Remarkable  are  the  strong  compressive  stresses  after 
heat  treatment  in  the  middle  (4  mm  depth)  both  for  the  z- 
direction  and  the  y-direction  (Fig.  4).  There  is  evidence 
that  solidification  partitioning  can  lead  to  a different  y/y'- 
microstructure  in  the  weld  zone,  thus  affecting  the  mean 
lattice  parameter  measured  by  diffraction.  The  immediate 
results  are  apparent  stresses  caused  by  changes  in  the 
unstressed  <f-spacing.  This  effect  can  be  removed  with 
some  confidence  only  for  the  measurements  close  to  the 


Fig.  4.  Residual  stresses  eryy  (parallel  to  the  weld)  before  and 
after  the  heat  treatment. 
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surface  ( = 1 mm  depth)  by  using  mechanical  equilibrium 
conditions.  Further  investigations  are  under  way  to 
clarify  this  effect  for  larger  depths. 

The  stresses  parallel  to  the  weld  (Fig.  4)  are  usually 
the  highest  and  the  most  dangerous  stresses  because  they 
reach  the  yield  limit  (around  600  MPa).  If  the  weld 
direction  is  oriented  such  that  the  weld  line  is  parallel  to 
the  main  stress  axis  - i.e.,  radially  for  a rotating  blade  - 
then  further  loading  would  exeed  the  How  stress  and 
cause  a dangerous  local  weakening  of  the  turbine  blade. 
Here,  the  effect  of  the  heat  treatment  is  most  beneficial 
because  almost  complete  stress  relief  is  achieved. 

The  stress  relief  effect  of  the  heat  treatment  extends 
also  to  the  microstresses  between  sub-grains  by  removing 
the  dislocations  created  by  the  plastic  deformation  in  the 
weld  region.  While  direct  observations  of  dislocations  are 


Fig.  5.  FWHM  at  the  best  instrument  resolution  before  and  after 
heat  treatment.  These  results  are  representative  for  all  measured 
depths. 


not  available,  we  can  use  the  full  width  at  half  maximum 
(FWHM)  of  the  measured  integrated  strain  distributions 
as  a measure  for  the  dislocation  density  (see  Fig.  5). 

Figure  5 shows  clearly  that  substantial  plastic  defor- 
mation exists  only  in  the  weld  itself  (0...  1 .0  mm)  and 
within  a very  narrow  zone  of  about  0.5  ...  1 .0  mm  in  the 
base  material.  This  effect  vanishes  completely  after  the 
heat  treatment. 

These  results  demonstrate  that  a heat  treatment  after 
welding  is  essential  for  avoiding  detrimental  effects  to 
the  mechanical  properties  of  the  turbine  blade. 


stresses:  y-direction 


EXPGUI:  Smoothing  the  Path  to  Powder 

Crystallography 





Back  in  the  dark  ages  of  computational  crystallogra- 
phy. cry  stallographers  were  expected  to  be  able  to  modify 
their  computer  programs  to  suit  the  computation  at  hand. 
Most  cry  stallographers  could  create  programs  to  convert 
data  between  different  formats  and  felt  comfortable 
reading  through  a program's  source  code  when  the 
program  documentation  did  not  adequately  explain  a 
topic.  (There  is  disagreement  on  when  this  ended,  some 
say  the  1950’s,  others  say  the  1 980 's  or  even  the  early 
2 1 st  century! ) 

The  advent  of  desktop  computing  has  also  introduced 
a new  paradigm  of  computing  — the  graphical  user 
interface  (GUI).  Visual  analogies  to  physical  devices 
(such  as  switches)  makes  operating  GUI  intuitive:  a good 
GUI  implementation  significantly  reduces  learning  time. 

Simultaneous  with  the  changes  in  computing,  crystal- 
lography has  evolved  from  a technique  used  exclusively 
by  specialists  to  one  of  many  tools  in  the  arsenal  of 
chemists,  physicists  and  engineers.  To  facilitate  this,  the 
NCNR  is  developing  GUI-based  software  for  crystallog- 
raphy. There  are  two  goals:  to  increase  the  utility  of  the 
softw  are  for  the  non-specialist  and  to  improve  the 
productivity  of  experts.  This  allow  s novices  to  concen- 
trate on  the  science  and  techniques  of  neutron  pow  der 
crystallography  rather  than  devote  their  efforts  to  learning 
arcane  details  of  computer  softw  are. 

The  Rictveld  refinement  technique  for  crystallo- 
graphic analysis  w as  adopted  within  the  NCNR  immedi- 
ately after  its  first  publication  ( 1 1 The  NCNR  then 
produced  Rictveld  software  specifically  adapted  to  the 
local  instrumentation,  which  is  still  being  improved  [2,3], 
However,  the  General  Structure  Analysis  System  (GSAS) 
package  of  programs  from  Los  Alamos  is  also  widely 
used  1 4 1 GSAS  has  wide  applicability  It  can  be  used 
w ith  virtually  any  ty  pe  of  neutron  or  x-ray  diffraction 
data,  including  both  powder  and  single  crystal  data.  It 
supports  a w ide  range  of  geometric  and  compositional 
restraints  and  constraints.  While  many  new  features  have 
been  added  to  GSAS  on  a continual  basis,  the  original 
text-based  user  interface  has  not  changed.  Novice  and 
occasional  users  often  find  the  many  levels  of  dialogs 
cumbersome  and  difficult  to  learn. 


Work  on  a GUI  for  GSAS  was  begun  several  years 
ago.  w ith  the  modest  goal  of  providing  a mechanism  for 
accessing  a small  portion  of  the  GSAS  features.  The 
project,  now  called  EXPGUI.  has  grown  in  scope  [5|. 
Most  of  the  commonly  used  pow  der  diffraction  features 
are  now  implemented,  allowing  a complete  Rictveld 
refinement  to  be  performed  graphically.  Features  not 
implemented  in  the  GUI  can  still  be  accessed  via  the 
standard  GSAS  user  interface.  EXPGUI  runs  identically 
on  all  platforms  where  GSAS  is  supported,  including 
both  Unix  and  Window  s In  FY2001.  the  Advanced 
Technology  Program  recognized  the  value  of  this  effort 
and  provided  funding  to  accelerate  development. 

The  GUI  is  designed  to  display  information  and 
controls  in  a concise  and  intuitive  fashion.  An  example 
screen,  shown  in  Fig.  1.  manipulates  structural  param- 
eters. This  figure  demonstrates  how  four  atoms  that  share 
sites  in  this  material  can  be  selected  together,  so  that  their 
refinement  options  can  be  changed  as  a group.  Another 
screen,  shown  in  Fig.  2,  applies  constraints  that  link 
appropriate  parameters  in  the  refinement. 
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Fig.  1.  Sample  screen  from  EXPGUI,  showing  where  atomic  pa- 
rameters are  displayed. 

Recent  EXPGUI  development  has  been  aimed  at 
extending  the  capabilities  of  the  GSAS  package.  As  an 
example,  w hen  EXPGUI  w as  recently  given  the  ability  to 
import  atomic  coordinates  from  files,  support  was  added 
for  CIF.  the  complex  standard  crystallographic  informa- 
tion file  format  developed  by  the  International  Union  of 
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Crystallographers.  Further,  the  import  routine  was  written 
in  a general  fashion,  so  that  other  formats  are  easy  to  add 

Where  possible.  EXPGUI  takes  advantage  of  the 
graphical  nature  of  modern  computers  to  offer  scientists 
better  access  to  refinement  tools.  Figure  3 shows  a 
graphical  display  of  the  observed  and  computed  diffrac- 
tion pattern.  Through  a simple  set  of  mouse  operations, 
the  user  has  "zoomed  in"  on  a portion  of  the  data,  to 
better  see  the  poor  agreement  between  observed  and 
computed  intensity  values,  indicating  that  the  structure  is 
not  well  modeled. 

This  graphical  tool  also  allows  the  user  to  see  the 
fitted  background  and  to  label  the  reflection  indices  — 
features  not  previously  available  within  GSAS.  Another 
example  is  a new  tool  for  background  fitting.  In  most 
cases,  background  fitting  is  nearly  automatic,  but  for 
some  refinements,  fitting  can  be  quite  difficult.  GSAS 
allows  a scientist  to  define  the  background  cither  through 
use  of  a function  with  rcfinable  terms  or  as  a spline 


Fig.  2.  EXPGUI  screen  for  entering  and  display  of  atomic  con- 
straints. In  this  case,  atoms  sharing  sites  are  constrained  to  have 
the  same  Ujso  values  and  have  fractional  occupancies  that  total 
to  unity. 

drawn  through  a set  of  fixed  points.  A new  tool  in 
EXPGUI  improves  GSAS  by  allowing  a user  to  define 
w here  the  background  curve  should  occur  using  a 
computer  mouse.  Since  many  experts  agree  that  fixed 
background  points  are  a poor  way  to  treat  background 
fitting.  EXPGUI  offers  an  augmented  feature:  a back- 
ground function  can  be  fit  to  the  user’s  desired  back- 
ground curve.  The  terms  for  this  function  will  be 


sufficient  to  obtain  an  initial  model  with  good  agreement 
to  the  data.  At  the  end.  the  background  terms  can  be 
optimized,  to  obtain  the  best  possible  model.  This 
capability  was  not  previously  av  ailable  in  GSAS. 

EXPGUI  has  been  widely  adopted  by  Rietveld  users. 
It  is  used  in  most  major  neutron  and  synchrotron  facili- 
ties. as  well  as  dozens,  if  not  hundreds,  of  universities, 
companies  and  research  centers  around  the  world. 
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Fig.  3.  Graphical  display  of  fit  in  EXPGUI.  This  plot  has  been 
“zoomed”  in  to  show  details.  Several  features  of  this  plot,  for 
example  the  display  of  reflection  indices  and  the  location  of  the 
fitted  background  were  not  previously  available  in  GSAS. 
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Neutron  Diffraction  Contributes  to  Improving  the 

Fatigue  Life  of  Bridges 


Welded  attachments  to  large-scale  structures  such  as 
bridges  can  limit  the  serviceability  and  prevent  effective 
use  of  higher  strength  steels.  To  avoid  this,  it  is  desirable 
to  enhance  the  fatigue  resistance  of  common  attachment 
details  such  as  transverse  stiffeners,  cover  plates,  gusset 
plates  and  other  welded  details  that  - at  the  weld  toe  - 
arc  often  the  initiation  point  for  crack  grow  th 


Fig.  1.  Partial  details  of  one  of  the  beam  specimens  being  fa- 
tigue tested  at  the  ATLSS  facility.  The  cover  plate  location  is 
shown  as  “W”  and  “Y”  on  the  right  side. 

Enhancement  of  fatigue  resistance  of  w elded  joints  by 
plastic  deformation  of  the  surface  and  by  improvement  of 
weld  toe  characteristics  is  well  established  1 1.2].  It  is 
known  that  the  conventional  improvement  techniques 
such  as  grinding,  shot  peening,  air  hammer  pcening.  gas 
tungsten  arc  (TIG)  re-melting,  and  use  of  improved 
welding  consumables  can  improve  fatigue  resistance  of 
welded  details  [3 1 . However,  these  procedures  arc 
manpower  intensive,  and  not  always  efficient  or  environ- 
mentally friendly.  Ultrasonic  Impact  treatment  (UIT) 
offers  an  alternative  means  to  avoid  the  negative  aspects 
of  other  methods  |4| 

The  post-w  eld  enhancement  of  w elded  details  by  UIT. 
developed  by  Statnikov  and  co-workers  [4|,  involves 
deformation  treatment  of  w elds  by  impacts  at  ultrasonic 
frequency  on  the  w eld  toe  surface.  The  objective  of  the 
treatment  is  to  introduce  beneficial  compressive  residual 
stresses  at  the  treated  weld  toe  zones  and  to  reduce  stress 
concentration  by  improving  the  weld  toe  profile.  UIT  is 
claimed  to  involve  a complex  effect  of  strain  hardening, 
reduction  in  welding  strain,  relaxation  in  residual  stress, 
and  a reduction  in  operating  stress  concentration,  thereby 


achieving  a deeper  cold  w orked  metal  I aver.  How  ever, 
virtually  none  of  the  effects  ascribed  to  UIT  have  been 
confirmed  by  appropriate  measurements. 

Recent  work  at  the  NCNR  has  made  significant 
contributions  in  providing  new  insights  into  how  UIT 
actually  affects  residual  stresses  in  the  weld  region 
around  a full-scale  cover  plate  for  a steel  girder.  In  Fig.  1 
a typical  I-beam  fatigue-tested  at  the  Center  for  Ad- 
vanced Technology  for  Large  Scale  Structures  (ATLSS) 
at  Lehigh  University  is  show  n.  In  the  initial  phase  of  the 
NCNR-Lehigh  collaboration,  neutron  diffraction  w as 
used  to  determine  the  effect  of  UIT  on  a simulated  flange 
plate  with  attached  cover  plate  detail  having  a 0.5  in. 
(12.7  mm)  end  fillet  weld  before  and  after  the  treatment 
by  UIT. 

The  geometry  and  dimensions  of  the  welded  detail  are 
the  same  as  in  the  full-scale  fatigue  test  beams  of  Fig.  1 
Both  base  plate  and  cover  plate  are  A572  Gr.50.  This 
kind  of  cover-plate  detail  is  largely  used  in  steel  bridges 
but  has  a low  level  of  fatigue  strength  in  comparison  w ith 
other  welded  configurations. 

The  neutron  diffraction  measurements  were  made  in  a 
fairly  complex  mesh,  mostly  in  the  near  vicinity  of  the 
weld  toe  as  shown  in  Fig.  2.  Overall,  residual  stresses 


Fig.  2.  Top  of  cover  plate  region  of  interest,  and  part  of  measure- 
ment mesh. 

were  determined  at  30  positions  in  the  weld  region, 
before  and  after  UIT  Gauge  volumes  of  3 x 3 x 3 mm3 
(interior)  and  2x2x2  mm3  (near  surface)  were  used. 

The  unstressed  cZ-spacing.  d w as  determined  from  a 
small  piece  of  the  starting  material.  The  depth  depen- 
dence of  X-direction  stresses  is  shown  in  Fig.  3.  The 
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near-surface  residual  stresses  determined,  nominally  I 
mm  below  the  actual  surface,  arc  shown  in  Fig.  4. 

To  our  knowledge  this  is  the  first  determination  of  the 
triaxial  stress  distribution  in  the  vicinity  of  a weld  treated 


Fig.  3.  X-direction  residual  stresses  vs.  depth  at  the  weld  center. 
Closed  symbols  are  pre-UIT.  Uncertainties  are  typically  ± 20  MPa. 
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with  UIT.  In  Fig.  4 both  the  anisotropy  of  stresses  around 
the  untreated  weld  metal,  and  the  additional  anisotropy 
introduced  by  UIT  are  seen  Before  UIT.  the  stresses  just 
beneath  the  weld  toe  arc  comparable  and  highly  tensile  in 
the  X and  Y directions,  which  reach  and  arc  beyond  the 
minimum  yield  stress  of  base  plate  material,  whereas 
those  normal  to  the  surface  are  generally  only  ~ 100  MPa 
tensile.  After  UIT  the  X-direction  stresses  become 
compressive,  changing  in  magnitude  by  about  400  MPa 
The  Y-direction  stresses  are  reduced  by  ~ 150  MPa  or 
less  but  remain  in  tension  except  at  one  point,  and  the 
normal  stresses  (Z-direction)  are  essentially  unchanged. 
The  UIT  effect  on  X-direction  stresses  (Fig.  3)  is  signifi- 
cant only  near  the  plate  surface,  but  deeper  than,  for 
example,  shot  peening.  The  depth  dependence  of  the  Y- 
and  Z-direction  stresses  w as  also  obtained 

Current  neutron  diffraction  studies  in  this  collabora- 
tion are  focusing  on  the  mapping  of  stress  gradients 
within  and  near  weld  joints  in  sections  of  actual  fatigue- 
tested  I-beams,  information  not  obtainable  by  any  other 
technique. 
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Fig.  4.  The  residual  stresses  determined  by  neutron  diffraction 
for  the  “near-surface”  points.  Uncertainties  due  to  counting  sta- 
tistics are  typically  about  ± 20  MPa.  Solid  symbols  are  before 
UIT. 
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Neutron  Imaging  and  Prompt  Gamma  Activation 
Analysis  Using  a Monolithic  Capillary  Neutron  Lens 


Glass  polycapillary  fibers  with  hundreds  of  hollow 
micrometer-wide  channels  have  been  used  for  guiding 
and  changing  the  direction  of  slow  neutrons  1 1 | The  use 
of  a monolithic  lens  for  PGAA  has  been  explored  earlier 
1 2 1 In  addition,  the  converging  nature  of  a focused  beam 
can  provide  magnified  images  for  neutron  radiography 
[3 1 We  have  combined  these  two  uses  of  the  focused 
beam  from  a monolithic  lens  and  have  implemented  an 
alignment  procedure  for  PGAA  aided  by  neutron  imag- 
ing. 

The  measurements  have  been  performed  at  the  end  of 
the  NG-0  beam  line  [4|  The  average  wavelength  of  the 
incident  beam  with  a Be  filter  is  about  0.6  urn.  Neutrons 
arc  imaged  w ith  a camera  that  consists  of  a ZnS(Ag)/6LiF 
scintillator  coupled  with  a CCD  camera  with  an  average 
pixel  size  of  50  pm.  The  gamma-ray  spectrometry  system 
[5 J consists  of  a p-type  germanium  photon  detector  w ith 
a virtual  pulser  "loss  free'  counting  module.  The  two 
monolithic  lenses  used  in  this  study  are  both  fabricated 
from  a bundle  of  capillaries  that  is  50  mm  long  and  w ith 
a hexagonal  cross  section  of  10  mm  w ide  ( flat-to-flat)  at 
the  entrance  and  tapered  to  3 mm  at  exit.  The  distance 
from  the  exit  to  the  focus  is  8.8  mm  for  lens  A and  9.8 
mm  for  lens  B Each  lens  is  installed  on  x-y-z  translation 
stages  w ith  two  directions  of  tilt  motion.  The  experimen- 
tal layout  is  shown  schematically  in  Fig.  I 


Z 

Neutron 

camera 


Fig.  1.  Schematic  experimental  layout. 

An  optimum  focal  spot  from  lens  A is  recorded  with 
the  imaging  detector.  Line  profiles  across  the  focal  spot 
in  both  the  horizontal  and  vertical  directions  are  taken  to 
determine  the  intensity  distribution  and  the  FWHM  (H: 


87  pm.  V:  100  pm)  of  the  focus  (Fig.  2 a).  The  gain  in 
neutron  current  density  within  the  100  pm  x 100  pm  area 
is  46  ± 2.  and  within  a 50  pm  x 50  pm  area  it  is  71  ± 5. 
Images  and  their  sizes  and  intensities  as  a function  of  z 
are  show  n in  Fig.  2b 
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Fig.  2.  (a)  Line  profiles  of  the  focal  spot  intensity  in  the  horizon- 
tal and  vertical  directions.  Insert:  Focused  beam  image  - the  focus 
is  near  the  center  of  the  hexagonal  shadow  cast  by  the  lens  en- 
trance, surrounded  by  the  background  of  unguided  neutrons, 
(b)  Intensity  and  base  width  of  the  spot  along  the  beam  path,  and 
images  of  the  focused  beam  “foot  print”.  A linear  fit  to  the  re- 
gion between  10  and  16  mm  gives  a slope  of  0.233  (rad),  or  a 
13.3°  convergence  angle. 

The  focused  beam  creates  a neutron  "point  source”, 
with  neutrons  departing  the  focus  within  a cone  of  about 
13°.  Any  neutron-absorbing  object  w ithin  the  path  of  the 
divergent  beam  w ill  cast  a magnified  shadow  at  sonic 
further  distance.  We  make  use  of  this  effect  in  order  to 
align  a small  sample,  a Gd  foil  (thickness  25  pm)  of  size 
2 mm  x 0.9  mm,  at  the  focus  for  PGAA  (Fig.  3a).  We 
vary  the  distance  between  the  sample  and  the  lens 
focus  by  translating  the  lens  along  the  direction  of  the 
beam.  The  magnified  image  enables  a better  alignment  of 
the  sample  with  respect  to  the  focus.  From  the  earlier 
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measurement  of  the  lens  focus  image,  we  determine  the 
location  of  the  focus  with  respect  to  the  large  hexagonal 
shadow  cast  by  the  lens  entrance.  We  then  move  the  lens 
such  that  the  image  of  the  sample  is  now  superimposed 
on  the  image  of  the  focus.  This  guarantees  that  the  focal 
point  is  indeed  on  the  sample.  Once  the  image  is  cen- 
tered, the  lens  is  translated  to  L ( ^ 0 such  that  the  focus  is 
on  the  sample,  ready  for  prompt  gamma  measurements. 


Fig.  3.  (a)  Magnified  image  of  the  sample  at  various  Lr  The  im- 
age of  the  sample  without  the  lens  is  also  shown  for  comparison, 
(b)  Prompt  gamma  spectrum  for  Gd  (182  keV)  and  for  In  (186 
keV)  measured  with  the  lens,  and  the  corresponding  image  of 
the  sample. 

At  the  center  of  the  same  Gd  sample  described  above, 
there  is  an  additional  indium  foil  of  about  I mm  \ 1 mm 
in  cross  section,  in  contact  with  the  Gd  strip.  The  prompt 
gamma  measurement  without  the  lens  gives  a peak 
intensity  for  Gd  at  182  keV  of  2.1  counts  per  second 
(cps),  but  the  In  peak  at  1 86  keV  is  not  visible  even  for  a 
20  min  measurement.  When  the  lens  is  used  in  a 10  min 
measurement,  not  only  has  the  Gd  peak  intensity  in- 
creased, but  the  In  peak  has  also  become  visible  (Fig. 

3b).  Since  the  Gd  piece  is  much  larger  than  the  focal  spot 
size,  the  gain  in  sensitivity  for  Gd  is  only  a factor  of  4. 
The  gain  for  In  is  not  estimated. 

We  have  also  used  lens  B for  an  even  smaller  sample 
a Gd  piece  of  128  pm  x 103  pm  x 25  pm  in  size  and 
2.59  pg  in  mass  (effective  mass  « 1 pg  taking  into 
account  self-absorption).  The  images  from  optical 
microscope  and  from  neutron  measurements  are  shown  in 
Fig.  4a.  Without  the  lens,  the  neutron  camera  cannot 
resolve  the  100  pm  particle.  The  gamma  intensity  with 
the  lens  is  17.4  cps  with  a 1.6  % la  counting  statistics 
after  1 hour  of  counting,  and  with  the  direct  beam  0.506 


cps  (3.9  % la)  after  19.5  hours.  The  gain,  given  by  the 
ratio  of  the  two  count  rates,  is  34  ± 1 .4.  A scan  of  the  lens 
in  the  vertical  direction  yields  a peak,  given  by  the 
convolution  of  the  sample  width  with  the  beam  width,  as 
shown  in  Fig.  4b.  The  FWHM  of  the  peak  is  148  ( ± 1.5) 
mm,  correctly  corresponding  to  the  100  pm  sample-width 
and  the  100  pm  beam-width  added  in  quadrature. 

The  effort  presented  here  is  a step  toward  providing 
quantitative  elemental  information  to  traditional  neutron 
radiography.  This  new  tool  will  provide  a unique  capabil- 
ity in  non-destructive  analysis  for  industrial  materials.  In 
the  future,  the  system  will  be  automated  to  scan,  giving  a 
more  efficient  probe  for  micro-analysis  with  neutron 
beams. 
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Fig.  4.  (a)  Optical  image  of  the  Gd  piece  in  comparison  with  a 
strand  of  human  hair  (left),  and  neutron  image  using  lens  B.  The 
large  solid  rectangle  is  a larger  Gd  strip  marker  to  help  locate 
the  smaller  sample  more  easily,  (b)  Gd  peak  intensity  produced 
by  a y-scan  of  the  lens. 
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Giant  Anharmonicity  and  Theory  of  Surprising  BCS 

Superconductivity  in  MgB2  at  40  K 


The  recent  surprise  discovery  of  superconductivity  in 
MgB,  at  39  k has  stimulated  a great  deal  of  research  on 
this  intercalated  graphite-like  system  (see  Fig.  1 ). 

Sparked  by  this  discovery,  we  set  out  to  unlock  the 
structural  secrets  and.  in  particular,  to  reveal  the  origin  of 
the  high  T in  MgB.:  an  electron-phonon  or  other  exotic 
mechanism?  To  answer  this  fundamental  question,  we 
calculated  T . its  pressure  dependence,  and  its  isotope 
effect  from  the  electronic  band  structure  and  lattice 
dynamics  of  MgB,  using  density  functional  theory  within 
the  generalized  gradient  approximation  [1-3]. 


Fig.  1.  The  crystal  structure  of  MgB2  consisting  of  B and  Mg  hex- 
agonal layers.  The  in-plane  boron  modes  (shown  by  arrows)  are 
strongly  coupled  to  the  boron  p>y  <r  bands  shown  as  the  green 
contour  and  isosurface  plots. 

Figure  2 shows  that  the  features  in  the  calculated 
phonon  density  of  states  (DOS)  are  in  excellent  agree- 
ment with  the  neutron  data  (GDOS).  giving  confidence 
that  the  calculations  provide  a sound  description  of  the 
physical  properties  of  the  system.  The  DOS  consists  of 
two  bands  of  phonons,  one  below  40  meV  corresponding 
primarily  to  acoustic  phonon  modes,  and  one  above  50 
meV  that  mostly  involves  the  boron  motions.  Inspection 
of  the  calculated  phonon-dispersion  curves  that  make  up 
the  high-energy  band  in  the  DOS  reveals  that  the  in-plane 
boron  phonons  (as  depicted  in  the  inset  to  Fig.  3)  are  very 
anharmonic.  To  demonstrate  this,  in  Fig.  3 we  plot  the 
total  energy  as  boron  atoms  move  according  to  one  of 
these  zone-center  in-plane  phonons  with  £ symmetry. 
The  potential  displays  a very  strong  anharmonic  term. 
Numerically  solving  the  Schrodinger  equation  for  this 
anharmonic  potential  yields  a phonon  energy  of  tico(E^) 


= 74.5  meV.  a 25  % enhancement  over  the  harmonic 
value  of  60.3  meV.  This  value  is  in  good  agreement  with 
recent  Raman  measurements.  The  giant  anharmonicity 
revealed  gives  the  first  hint  that  the  in-plane  modes  are 
strongly  coupled  to  the  p a bonding  orbitals  of  boron, 
as  shown  schematically  in  Fig.  I . This  coupling  is  also 
evident  from  the  splitting  of  the  boron  a bands  (red  lines) 


Fig.  2.  Generalized  (top)  and  the  calculated  (bottom)  phonon 
density  of  states.  The  intensities  of  the  peaks  do  not  agree  well 
because  in  DOS  B and  Mg  ions  contribute  equally  while  in  GDOS 
they  are  weighted  by  neutron  cross-sections  and  inverse  masses. 

with  the  E,r  phonons  (see  Fig.  3).  Note  that  the  other 
bands  are  not  affected  by  the  E1  phonons. 

The  splitting  of  the  boron  a bands,  when  averaged 
over  the  Fermi  surface,  gives  an  electron-phonon  (EP) 
coupling  constant  A « 1 . Using  this  value  in  the  McMillan 
expression  for  T with  m (£v)  and  taking  a typical  value 
for  the  Coulomb  repulsion  //*  = 0. 1 5,  we  obtain  a T of 
39.4  k,  in  excellent  agreement  with  experiments.  We  also 
solved  the  Schrodinger  equation  for  the  potential  shown 
in  Fig.  3 for  different  boron  masses  and  obtained  «AE,)  = 
291.8  M'0515 and  A = 0.615!  A/0169,  which  yields  the  T -M 
curve  shown  in  Fig.  4 and  a boron  isotope  effect  a = 

0.21.  in  good  agreement  with  the  experimental  value  of 
0.26  ± 0.03. 

Since  the  pressure  dependence  of  T puts  a stringent 
test  on  any  theory  of  superconductivity,  we  repeated  the 
calculations  of  phonons  and  electronic  band  structures 
discussed  above  for  isotropic,  uniaxial  (along  c-axis),  and 
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biaxial  (in  the  a^-plane)  pressures  [3],  We  find  that  while 
co(Ej  increases  with  increasing  pressure,  the  density  of 
states  at  the  Fermi  energy  decreases.  The  EP  constant  A 
also  shows  significant  changes  with  pressure.  Inserting 
all  these  competitive  effects  into  the  McMillan  formula 
yields  the  pressure  dependence  of  T shown  in  Fig.  4. 

For  isotropic  pressure.  T decreases  with  increasing 
pressure  almost  linearly  at  a rate  of  « -1.0  K/GPa,  in 
excellent  agreement  with  the  experimental  value  of -1.1 
K/GPa.  We  also  predict  a cusp  in  the  T -P  curve  around  P 
« 20  GPa.  A similar  cusp  was  recently  observed  experi- 
mentally. Our  calculations  indicate  that  7 should  in- 


Boron  displacement  in  E Frozen  Phonon  (A) 


Fig.  3.  Top:  Total  energy  curve  as  a function  of  boron  displace- 
ment for  E2g  mode  (shown  in  the  inset),  indicating  a large 
anharmonic  term  in  the  potential.  Bottom:  Band  structure  of  the 
undistorted  (left)  and  distorted  structures  (right)  by  E2g  phonons 
(us  * 0.06  A).  See  Ref.  [1]  for  the  animation  of  the  zone  center 
phonons  and  their  coupling  with  the  B a bands. 


crease  first  and  then  decrease  with  increasing  ocompres- 
sion,  while  it  should  decrease  rapidly  with  c/A-compres- 
sion.  Hence,  when  single  crystal  samples  of  MgB, 


Fig  4.  Top:  Boron  mass  dependence  of  the  T.  Bottom:  Pressure 
dependence  of  T as  a function  of  uniform,  ab- axes,  and  c-axis 
compression,  respectively. 

become  available,  measurements  of  the  ab-  and  c- 
compression  dependence  of  T should  provide  a critical 
test  of  our  theory. 

MgB.  may  be  the  ultimate  BCS  s-wave  superconduc- 
tor. with  parameters  controlling  T fully  optimized  to 
yield  the  highest  possible  T . However,  even  if  T cannot 
be  increased  further,  the  low  cost,  light  mass,  easy 
fabrication,  nearly  isotropic  high  conductivity  of  MgB., 
which  also  has  a large  critical  current,  will  no  doubt  find 
many  important  technological  applications  in  the  near 
future. 
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Polaron  Formation  and  Colossal 

Magnetoresistance  in  Manganites 


The  recent  discovery  of  huge  magnetoresistance 
effects  in  the  manganese  oxide  class  of  materials  (such  as 
La:  A MnO,  (A=Sr,  Ca,  Ba) ) has  rekindled  intense 
interest  in  these  systems  for  three  reasons.  One  is  that 
this  “colossal  magnetoresistance”  (CMR)  offers  potential 
in  a number  of  technologies,  such  as  for  read/write  heads 
in  magnetic  recording  media,  sensors,  and  spin-polarized 
electronics.  The  second  reason  is  that  these  materials  are 
interesting  from  a fundamental  viewpoint.  In  contrast  to 
traditional  isotropic  ferromagnets  such  as  Fe,  Co,  Ni,  and 
EuO.  where  the  spin  system  is  isolated  from  the  lattice,  in 
the  manganites  the  charge,  spin,  and  lattice  degrees  of 
freedom  are  strongly  coupled  together.  This  leads  to  a 
delicate  balance  of  interactions  that  gives  rise  to  a rich 
variety  of  physical  phenomena  of  current  interest  in 
condensed  matter  science.  These  include  a metal- 
insulator  transition  that  is  juxtaposed  with  ferromag- 
netism, charge  ordering,  orbital  ordering,  polaron 
formation,  electronic  phase  separation,  and  spin  and 
charge  stripes.  The  third  reason  is  because  of  the  similar- 
ity to  the  high  T cuprate  oxides.  Recent  progress  in  our 
understanding  of  the  cuprates  has  provided  insights  into 
the  manganites,  and  a deeper  understanding  of  the 
fundamental  properties  of  the  manganites  will  surely 
elucidate  the  shared  concepts  underlying  both  classes  of 
materials. 

The  electrical  conductivity  in  these  manganites  is 
directly  connected  to  the  magnetic  system  through  the 
double  exchange  mechanism,  while  the  Jahn-Teller 
distortion  of  the  Mn,+Oh  couples  the  magnetic  and  lattice 
systems,  resulting  in  the  electronic,  lattice,  and  magnetic 
degrees  of  freedom  being  intimately  intertwined.  The 
ferromagnetic  transition  then  is  non-conventional,  and 
typically  coincides  with  a metal-insulator  transition  to 
produce  CMR.  The  present  measurements  on 
La07Can  ,MnO,,  which  is  near  optimal  doping  for  CMR, 
show  that  polarons  form  as  the  ferro-  to  paramagnetic 
transition  is  traversed,  self-trapping  the  carriers.  These 
polarons  then  “melt”  when  a magnetic  field  is  applied, 
freeing  the  carriers  and  generating  the  CMR. 

The  scattering  from  a lattice  polaron  arises  from  the 
structural  distortion  that  surrounds  a carrier  and  traps  it. 


Fig.  1.  Polaron  correlation  peaks  near  the  5 0 0 Bragg  peak  at 
270  K,  just  above  Tc  = 255  K. 

Individual  polarons  generate  diffuse  scattering  around  the 
fundamental  Bragg  peaks,  while  polaron  correlations  will 
give  rise  to  additional  superlattice  peaks.  We  observe 
both  of  these  effects  [1],  and  in  Fig.  1 we  show  the  four 
satellite  peaks  at  'A  1 A 0 positions  near  the  5 0 0 Bragg 
peak,  just  above  T (255  K).  The  intrinsic  widths  of  these 
peaks  indicate  that  the  correlations  are  short  range  («  15 
A)  in  nature,  and  this  range  is  found  to  be  only  weakly 
temperature  dependent.  Polarized  neutron  measurements 
show  that  these  peaks  originate  from  lattice  distortions 
rather  than  from  magnetic  scattering.  Similar  direct 
evidence  for  polaron  formation  and  correlations  has  been 
observed  in  the  bilayer  manganite  Lat  ,Sr  gMn,07  [2], 
although  the  ordering  wave  vectors  and  the  ratios  of 
diffuse  scattering  to  correlated  scattering  are  different. 

A model  that  can  explain  both  the  ordering  wave 
vector  as  well  as  the  overall  behavior  of  the  observed 
intensities  of  the  short-range  correlation  peaks  is  the  CE 
model  developed  for  the  half-doped  case.  The  CE  model 
has  an  orbital ly-ordered  Jahn-Teller  lattice  with  charge 
stripes  in  the  [110]  direction,  and  an  antiferromagnetic 
ground  state.  Such  ordering  has  been  observed  for  x > 0.5 
at  low  temperatures.  This  behavior  contrasts  with  the 
present  ferromagnetic  system,  where  the  orbital/charge 
correlations  are  short  range  in  nature,  occur  only  above 
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Tc,  and  are  not  accompanied  by  any  antiferromagnetic 
correlations.  The  CE  “stripes”  are  also  consistent  with 
stripe  formation  as  seen  in  other  transition  metal  oxides 
(nickelates  and  cuprates). 

Of  particular  interest  is  how  the  polaron  ordering  is 
related  to  the  metal-insulator  transition  and  the 
quasielastic  “central  peak”  [3|  associated  with  the  spin 


Fig.  2.  Temperature  dependence  of  the  intensity  of  the  polaron 
peak  at  3.75  0.25  0,  compared  to  the  magnetic  quasielastic  cen- 
tral peak,  and  the  resistivity.  The  data  have  been  scaled  vertically 
so  the  peak  heights  match.  The  similarity  of  the  data  indicates  a 
common  physical  origin. 

component  of  the  polaron.  As  we  see  in  Fig.  2.  the 
temperature  dependence  of  the  intensity  of  the  charge 
ordering  is  nearly  identical  to  the  quasielastic  magnetic 
peak,  and  both  follow  the  resistivity  behavior.  This  is  a 
very  strong  indication  that  these  properties  have  a 
common  origin  in  polaron  formation. 

A further  indication  of  the  strong  relationship  between 
charge  ordering  and  resistivity  is  established  by  measure- 
ments in  magnetic  fields  [1],  In  Fig.  3 we  see  that  the 
intensity  of  the  charge  ordering  peaks  decreases  rapidly 
with  increasing  field.  Thus  the  polarons  “melt”  when  the 
magnetization  develops  and  the  conductivity  increases, 
either  by  lowering  the  temperature  (Fig.  2)  or  increasing 
the  applied  field,  producing  CMR.  Identical  field  depen- 
dence has  been  observed  both  for  the  central  peak  in 


LaQ7Can ^MnO,  [1,3]  as  well  as  the  polaron  peaks  in  the 
bi layer  system  [2]. 

The  above  experimental  results  reveal  that  both  the 
spin  and  charge  correlations  associated  with  the  polarons 
in  La()7Ca01MnO,  appear  together,  and  have  a very  similar 
temperature  dependence.  This  coincidence  may  explain 
the  amplified  magnetoresistive  effects,  as  well  as  the 
absence  of  conventional  magnetic  critical  behavior.  The 
physical  quantities  that  control  the  characteristic  tem- 
perature for  polaron  formation,  including  the  role  that 
intrinsic  inhomogeneities,  stress,  and  phase  separation 
play  in  their  formation,  still  remain  to  be  elucidated. 
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Fig.  3.  Magnetic  field  dependence  of  the  intensity  of  the  polaron 
peak,  showing  that  the  polarons  melt  with  increasing  field. 
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Ferromagnetism  and  Spontaneous  Vortex 
Formation  in  Superconducting  ErNi2B2C 


The  competition  between  magnetic  order  and  super- 
conductivity has  had  an  interesting  history.  The  first 
examples  of  true  long  range  magnetic  order  coexisting 
with  superconductivity  were  provided  by  the  ternary 
Chevrel-phase  superconductors  and  related  compounds. 
The  magnetic  ordering  temperatures  were  all  low  («  I K), 
indicating  that  dipolar  interactions  dominate  the  energet- 
ics of  the  magnetic  system.  Similar  behavior  is  found  for 
the  rare  earth  ordering  in  the  hole-doped  cuprates  [1  j.  The 
first  systems  where  exchange  clearly  dominates  the 
magnetic  energetics  were  provided  by  the  electron-doped 
cuprates  [ 1 ],  followed  by  the  borocarbides  of  central 
interest  here  [2], 

Among  these  systems,  the  rare  and  interesting  situa- 
tion where  ferromagnetic  interactions  are  present  has 
attracted  special  attention  because  of  the  competitive 
nature  between  the  superconducting  screening  and  the 
internally  generated  magnetic  field.  In  the  Chevrel  phase 
materials  this  competition  gives  rise  to  long  wavelength 
oscillatory  magnetic  states  and/or  reentrant  superconduc- 
tivity. Recently,  superconducting  ErNi,B,C  {T  = 11  K) 
appeared  to  develop  a net  magnetization  below  2.3  K, 
well  below  the  onset  of  long  range  spin-density-wave 
order  at  Tn=  6 K.  We  have  carried  out  a comprehensive 
study  of  the  magnetic  order  and  vortex  structure  in  this 
material,  and  directly  demonstrated  that  this  transition 
does  indeed  correspond  to  the  development  of  a net 
atomic  magnetization  that  coexists  with  superconductiv- 
ity, and  this  results  in  the  spontaneous  formation  of 
vortices  in  the  system  [3]. 

Superconducting  ErNi,B,C  orders  magnetically  at  6 K 
into  a transversely  polarized  spin  density  wave  structure 
[2,3],  with  the  modulation  wave  vector  d along  the  a axis 
and  the  spins  along  the  b direction.  A portion  of  a scan 
along  the  a axis  is  shown  in  Fig.  1 at  three  temperatures. 
The  initial  ordering  is  a simple  spin  density  wave  as 
shown  in  the  schematic,  with  an  incommensurate  modula- 
tion of  5 s 0.55  a * (a*  = 2n/a).  At  2.4  K additional 
higher-order  peaks  become  observable  as  the  magnetic 
structure  squares  up,  as  indicated  schematically  at  the 
bottom  of  Fig.  I . These  peaks  are  all  odd-order  harmonics 
as  expected  for  a square-wave  magnetic  structure. 


Below  the  weak  ferromagnetic  transition  (7WFM)  at  2.3 
K.,  the  data  (Fig.  1 ) indicate  that  a new  series  of  peaks  has 
developed,  which  are  even  harmonics  of  the  fundamental 
wave  vector.  Polarized  neutron  measurements  unambigu- 
ously establish  that  both  the  odd-order  peaks  and  the 
even-order  peaks  are  magnetic  in  origin. 

The  integrated  intensities  of  the  odd-  and  even-order 
harmonics  are  shown  in  Fig.  2 as  a function  of  T.  The 
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Fig.  1.  Scan  along  the  [hOO]  direction  at  7"  = 1.3  K,  2.4  K,  and 
4.58  K.  The  data  have  been  offset  for  clarity.  Above  the  weak 
ferromagnetic  transition  at  2.3  K we  observe  the  fundamental 
peak  at  ^ 0.55  along  with  odd-order  harmonics,  while  below  a 
new  set  of  even  harmonics  develops  (arrows).  Also  shown  is  a 
schematic  of  the  spin  density  wave  order  for  T just  below  7"N, 
and  the  (undistorted)  wave  when  it  squares  up. 

even-order  peaks  abruptly  develop  below  2.3  K,  con- 
comitant with  the  development  of  the  new  (weak) 
ferromagnetism  at  There  is  a substantial  thermal 
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hysteresis  associated  with  the  weak  ferromagnetic 
transition,  suggesting  that  this  transition  is  first  order  in 
nature. 

The  present  neutron  results  demonstrate  that  a net 
magnetization  develops  in  ErNi,B,C  in  the  magnetically 
ordered  state  at  low  temperatures,  making  this  the  first 
such  “ferromagnetic-superconductor”  since  HoMo6Sg, 
HoMo^Se^  and  ErRh4B4.  For  ErNi,B,C  the  net  magneti- 
zation is  much  smaller  than  for  these  earlier  systems, 
which  allows  coexistence  with  superconductivity  over  an 
extended  T range.  This  presents  the  intriguing  possibility 
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Fig.  2.  Integrated  intensities  of  the  odd  (open  circles)  and  even 
(solid  circles)  order  harmonics  as  a function  of  T. 

that  in  an  applied  field  vortices  will  form  spontaneously 
when  this  net  atomic  magnetization  is  present.  The  small 
angle  scattering  data  of  the  vortex  structure  (Fig.  3 top) 


Fig.  3.  Intensity  of  the  vortex  scattering  vs.  wave  vector  Q at 
850  Oe  above  and  below  the  weak  ferromagnetic  transition.  The 
shift  in  the  peak  position  demonstrates  that  additional  vortices 
spontaneously  form  at  low  temperatures.  The  temperature  de- 
pendence (bottom)  shows  that  this  spontaneous  vortex  formation 
is  directly  related  to  the  weak  ferromagnetic  transition. 

phases,  which  may  prove  useful  in  high  current  applica- 
tions of  superconductors. 


show  that  the  lattice  has  the  expected  spacing  at  lower 
applied  fields.  At  higher  fields,  as  the  field  begins  to 
penetrate,  vortices  spontaneously  form  in  addition  to 
those  expected  from  the  applied  field  alone,  increasing 
the  vortex  density  and  shifting  the  vortex  peak.  The  T 
dependence  of  this  shift  (Fig.  3 bottom)  makes  it  clear 
that  this  behavior  is  directly  related  to  onset  of  the  net 
magnetization  in  the  system.  It  is  also  interesting  to  note 
that  the  vortex  pinning  is  enhanced  in  both  magnetic 
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Structure  Of  A Local  Spin  Resonance  In  A 
Geometrically  Frustrated  Antiferromagnet 


For  many  crystalline  magnetic  materials,  only  a long 
range  ordered  spin  configuration  can  satisfy  all  near 
neighbor  spin  interactions.  Such  systems  generally 
display  a finite  temperature  transition  to  a broken 
symmetry  phase  with  long-range  magnetic  order.  How- 
ever, on  certain  lattices  with  weak  connectivity  and  a 
triangular  motif,  short-range  interactions  can  be  satisfied 
without  long-range  order  [1], 

To  explore  this  possibility  we  examined  magnetic 
order  and  fluctuations  in  ZnCro04.  The  B-site  of  this 
spinel  lattice  is  solely  occupied  by  spin-3/2  Cr3+  and  this 
leads  to  a magnet  with  dominant  nearest  neighbor 
interactions  on  the  lattice  of  corner-sharing  tetrahedra 
shown  in  Fig.  1 . 


Fig.  1.  Lattice  of  corner-sharing  tetrahedra  formed  by  the  octa- 
hedrally  coordinated  B sites  of  the  spinel  AB204  structure. 


Exchange  interactions  in  this  structure  are  satisfied  at 
the  mean  field  classical  level  when  every  tetrahedron 
carries  zero  net  magnetization.  Not  only  are  there  many 
ways  that  a tetrahedron  can  have  no  moment,  but  there 
are  also  many  ways  to  organize  such  non-magnetic 
tetrahedra  on  the  corner-sharing  lattice.  The  result  is  a 
“geometrically  frustrated"  spin  system  with  many  ways  to 
satisfy  interactions  without  long-range  order. 

Analysis  of  specific  spin  models  on  the  B-site  spinel 
(or  pyrochlore)  lattice  indicates  that  spin- 1/2  and  spin-oo 


models  have  short-range  order  down  to  T=  0 [2],  while 
long-range  order  is  induced  by  quantum  fluctuations  for 
an  unknown  intermediate  range  of  spin  values.  Experi- 
ments indicate  that  spin-3/2  ZnCr,04  is  “close”  to  the 
quantum  critical  point  that  separates  the  low  spin  quan- 
tum-disordered phase  from  the  intermediate  spin  long- 
range  ordered  phase.  Specifically,  the  relaxation  rate  for 
magnetic  excitations,  A follows  a power-law  that  ex- 
trapolates to  zero  as  T— > 0 indicating  quantum  criticality 
[3],  This  state  of  affairs  however  does  not  persist  to  the 
lowest  temperatures.  Instead,  a first  order  transition  at  Ts 
= 12.5  K causes  a transition  from  the  cubic  cooperative 
paramagnet  to  tetragonal  Neel  order. 


Fig.  2.  Constant-Q  scans  above  and  below  fN. 

Figure  2 shows  that  a gapless  continuum  of  magnetic 
scattering  above  TN  is  pushed  into  a local  spin  resonance 
at  tico~  4.5  meV  « ./with  remarkably  little  dispersion 
throughout  the  Brillouin  zone.  The  result  is  analogous  to 
the  spin-Peierls  transition  of  the  uniform  spin- 1/2  chain, 
where  quantum  critical  fluctuations  are  pushed  into  a 
finite  energy  singlet-triplet  transition  through  structural 
dimerization.  Our  recent  synchrotron  X-ray  and  neutron 
powder  experiments  indicate  that  deformation  of  tetrahe- 
dra does  indeed  occur  for  ZnCi\04.  However,  the  struc- 
tural changes  push  the  system  to  long  range  order  rather 
than  to  quantum  disorder,  as  indicated  by  the  magnetic 
Bragg  peaks  and  spin  wave  excitations. 
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Our  single  crystal  experiment  [4]  also  enabled  unique 
insight  into  the  local  structure  of  spin-fluctuations  in 
geometrically  frustrated  systems.  Figure  3 shows  the  Q- 
dependence  of  low  energy  magnetic  scattering  in  two 
high  symmetry  planes  above  and  below  Ts.  While  the 
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spectrum  for  spin  fluctuations  changes  dramatically  at  the 
first  order  phase  transition,  the  structure  factor  clearly 
does  not.  Also  shown  in  the  figure  is  the  structure  factor 
for  six  spins  of  <1 1 1>  type  kagome  hexagons  precessing 
with  n phase  shift  between  neighbors.  The  proposal  by  O. 
Tchernyshyov  et  al.  [5]  that  these  are  the  dominant  low 
energy  spin  fluctuation  for  spins  on  the  B-site  spinel 
lattice  is  clearly  borne  out  by  the  data. 

The  present  data  for  ZnCi\04  show  that  geometrically 
frustrated  lattices  have  composite  low  energy  degrees  of 
freedom  analogous  to  rigid  unit  modes  in  certain  open 
framework  lattice  structures.  To  better  understand  the 
unusual  type  of  phase  transition  that  occurs  in  this 
system,  it  must  be  determined  what  defines  the  4.5  meV 
energy  scale  for  hexagon  excitation  in  the  ordered  phase. 
Do  quantum  fluctuations  play  a significant  role  or  does 
the  broken  symmetry  between  exchange  interactions 
within  tetrahedra  induce  the  resonance?  The  answer  to 
this  question  is  now  being  pursued  through  an  accurate 
determination  of  the  complex  low  temperature  lattice  and 
magnetic  structure. 
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Fig.  3.  Color  images  of  inelastic  neutron  scattering  intensities 
in  the  (hkO)  and  (hkk)  symmetry  planes.  From  top  to  bottom  are 
shown  data  from  below  and  above  TN,  and  the  structure  factor 
for  the  hexagon  mode  proposed  by  O.  Tchernyshyov.  T = 1.7  K 
data  were  taken  integrating  over  the  energy  range  Tioj  from  3.0 
meV  to  6.7  meV.  Data  in  the  paramagnetic  phase  (T  = 15  K)  were 
taken  at  ftco  = 1 meV. 
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Pressure-Induced  Phase  Transition  of  C E 

Micelles 


j-  . | - : I — 

Nonionic  surfactants  form  a variety  of  microstructures 
in  water,  ranging  from  simple  micelles  at  low  surfactant 
concentrations  to  complex  mesophases,  such  as  hexago- 
nal or  lamellar  phases  at  high  concentrations.  The  effect 
of  pressure  on  the  structure  of  micelles  and 
microemulsions  has  not  been  extensively  studied.  Small 
angle  neutron  scattering  (SANS)  is  particularly  well 
suited  for  such  measurements  because  the  range  of  length 
scales  probed  includes  both  the  particle  size  and  the 
interparticle  spacing.  Pressure  effects  observed  in  SANS 
measurements  of  surfactant  microstructure  are  typically 
interpreted  in  terms  of  geometric  packing  arguments  that 
focus  on  the  compressible  hydrophobic  tails  of  the 
surfactants.  In  this  context,  increasing  pressure  has  the 
single  effect  of  decreasing  the  surfactant  tail  volume, 
thereby  increasing  the  curvature  of  the  oil-water  inter- 
face. Conversely,  increasing  temperature  dehydrates  the 
nonionic  headgroups,  decreasing  the  headgroup  area  and 
decreasing  the  film  curvature.  Thus,  temperature  and 
pressure  can  be  viewed  as  thermodynamic  variables  with 
inherently  different  mechanisms  and  opposite  effects  for 
controlling  microstructure  in  nonionic  surfactant  solu- 
tions. 

Here  we  report  the  results  of  high-pressure  SANS 
experiments  on  a solution  of  1 % by  weight 
pentaethylene  glycol  mono-n-dodecyl  ether  (CpE.)  in 
D,0  at  20  °C  and  pressures  up  to  « 300  MPa.  The  phase 
diagram  for  this  system  [ 1 ] at  ambient  pressure  is  shown 
in  Fig.  1 . At  this  temperature  and  surfactant  concentration, 
a single-phase  micellar  solution  (L,  phase)  forms  at 
ambient  pressure,  well  below  the  lower  critical  solution 
temperature  (LCST)  for  liquid-liquid  equilibrium  and  far 
removed  from  the  H:  hexagonal  phase  at  much  higher 
CpE5  concentrations.  The  microstructure  of  the  L^phase 
at  ambient  pressure  is  known  to  be  a network  of  branched 
semi-flexible,  cylindrical  micelles  with  the  branch  points 
comprised  of  three-armed  junctions. 

SANS  spectra  were  measured  using  the  NIST  high- 
pressure  cell  and  neutrons  of  wavelength  X = 6 A, 
covering  a c/-range  of  0.0 1 2 A'1  < q < 0.22  A'1 . The 
scattering  curves  obtained  at  3.4,  241,  and  255  MPa  are 
shown  in  Fig.  2.  The  curves  at  3.4  MPa  and  241  MPa  are 


virtually  identical,  indicating  no  significant  change  in 
microstructure  with  increasing  pressure  up  to  241  MPa. 
Fitting  these  curves  using  a form  factor  for  cylindrical 
micelles  gives  a radius  of  (2 1 .0  ± 0.2)  A and  a length 
greater  than  600  A.  independent  of  pressure.  However,  a 
small  increase  in  pressure  from  241  MPa  to  255  MPa 
leads  to  the  appearance  of  a peak  in  the  scattering 
intensity  at  q ~ 0.130  A-1,  indicative  of  a locally  ordered 
system.  A similar  transition  has  been  reported  in  high- 
pressure  SANS  studies  of  tetradecyldimethylaminoxide 
(TDMAO)  micelles  in  D,0  at  pressures  up  to  300  MPa, 
but  the  high-pressure  microstructure  was  never  deter- 
mined [2]. 


Fig.  1.  Temperature-composition  phase  diagram  for  C12E5/water 
solutions  at  ambient  pressures  [1]. 

Shown  in  Fig.  2 is  the  scattering  curve  for  the  H| 
hexagonal  phase  at  49  % by  weight  CpF5  at  20  °C,  and 
ambient  pressure.  The  peak  at  q « 0.120  A'1,  arising  from 
the  hexagonal  lattice  of  cylindrical  micelles,  is  similar  to 
the  peak  for  the  1 % by  weight  CpE5  solution  at  255 
MPa,  suggesting  that  this  new  high-pressure  phase  may 
resemble  a slightly  compressed  state  of  the  llt  hexagonal 
phase  at  ambient  pressure.  The  formation  of  a lamellar 
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phase  from  cylindrical  micelles  is  unlikely,  since  this 
corresponds  to  increasing  the  hydrophobic  core  volume- 
to-surface-area  ratio  per  surfactant  molecule,  or  equiva- 
lently, decreasing  the  spontaneous  curvature  of  the 
surfactant  film.  The  application  of  pressure  would  have 
the  opposite  effect.  We  conclude,  therefore,  that  the 
observed  change  in  microstructure  corresponds  to  a 
pressure-induced  L^-H,  phase  transition  from  a network 


Fig.  2.  Measured  SANS  spectra  at  20  °C  for  1 % by  weight  C12E5  in 
D20  as  a function  of  pressure  and  49  % by  weight  C12E5  in  D20  at 
ambient  pressure. 

of  branched  semi-flexible,  cylindrical  micelles  to  hexago- 
nal ly  ordered  bundles  of  cylindrical  micelles. 

To  further  understand  the  SANS  results,  we  have 
measured  the  temperature  dependence  of  the  L1-H| 
transition  pressure  and  find  that  the p-T curve  follows  the 
p-T  freezing  curves  for  liquid  rt-alkanes  of  comparable 
hydrocarbon  chain  length.  TWdecane  solidifies  at  a 
pressure  of « 250  MPa  at  20  °C,  which  is  close  to  the 
pressure  for  the  observed  1^-^  phase  transition.  We 
propose  that  the  CpE5  micelle  hydrophobic  core,  equiva- 
lent to  u-decane,  does  solidify  at  these  conditions,  such 
that  the  micelles  lose  flexibility,  and  hence  conforma- 
tional entropy.  An  analysis  of  the  geometric  packing 


constraints  for  three-arm  junctions  coexisting  with 
cylinders  shows  that  when  the  surfactant  tail  volume 
decreases  with  increasing  pressure,  the  fraction  of 
surfactant  forming  junction  points  also  decreases.  This 
indicates  that  the  formation  of  three-arm  junctions 
becomes  increasingly  unfavorable  at  higher  pressures  due 
to  the  compression  of  the  CpE.  micelle  hydrophobic 
core.  Our  calculations  predict  that  no  junctions  should  be 
present  at  P > 275  MPa  at  20  °C,  which  is  in  good 
agreement  with  our  observation  of  a structural  transition 
between  241  MPa  and  255  MPa  at  this  temperature. 
Consequently,  the  network  of  branched  threadlike 
micelles  becomes  globally  unstable.  The  formation  of 
hexagonally  ordered  bundles  of  cylindrical  micelles 
follows  as  the  attractive  van  der  Waals  forces  between 
the  micelles  are  not  offset  by  entropic  repulsive  undula- 
tion interactions  that  are  not  present  in  the  now-solidified 
hydrophobic  core. 

The  practical  significance  of  these  results  is  to  show 
that  pressure  allows  access  to  regions  of  the  CpE./water 
phase  diagram  that  are  virtually  inaccessible  to  tempera- 
ture. Thus,  the  use  of  pressure  may  offer  unique  ap- 
proaches for  directing  and  stabilizing  certain  surfactant 
microstructures  that,  in  turn,  could  prove  useful  for 
creating  novel  soft  materials. 
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Probing  The  Structure  Of  Aerosol  Nanodroplets 


Nanodroplet  aerosols  form  readily  in  the  supersonic 
expansions  that  occur,  for  example,  in  turbomachinery, 
jet  exhausts,  and  volcanic  eruptions.  Thus,  understanding 
particle  formation  and  growth  when  cooling  rates 
approach  106  K/s  is  of  broad  scientific  interest. 

From  a fundamental  point  of  view,  particles  with  radii 
< 10  nm  are  important  because  they  lie  in  the  critical 
transition  zone  between  large  molecular  clusters  and  bulk 
material.  In  multicomponent  droplets,  there  is  also  strong 
theoretical  evidence  for  surface  enrichment,  i.e.,  that  the 
surface  and  interior  compositions  can  differ  significantly. 
Surface  enrichment  is  important  because  it  affects 
nucleation,  growth  and  evaporation  kinetics,  and  the 
heterogeneous  chemistry  of  aerosol  droplets.  A major 
goal  of  our  work  is  to  use  small  angle  neutron  scattering 
to  find  direct  evidence  for  surface  enrichment  in 
nanodroplets. 


Fig.  1 . The  nozzle  is  placed  in  the  sample  chamber  at  right  angles 
to  the  neutron  beam. 

We  produce  nanodroplets  with  radii  between  5 nm 
and  12  nm  and  number  densities  in  the  range  ot  10"  cm  ’ 
to  1 O'2  cm2  by  rapidly  expanding  a dilute  vapor  mixture 
of  D,0  (or  other  condensible)  and  N,  in  a supersonic 


nozzle  apparatus  (Fig.  1 ) [1],  The  volume  fraction  of 
droplets  is  always  less  than  10'5  and,  thus,  the  scattering 
signals  are  close  to  the  instrument  background.  Further- 
more, the  high  velocity  of  the  droplets  (400  m/s  to 
500  m/s)  is  comparable  to  the  speed  of  the  neutrons.  This 
leads  to  a Doppler  shift  in  the  momentum  of  the  scattered 
neutrons  [2,3],  and  the  2-D  scattering  patterns  are 
anisotropic  (Fig.  2)  even  though  spherical  droplets  at  rest 
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Fig.  2.  The  observed  2-D  scattering  pattern  for  an  aerosol  with 
an  average  particle  velocity  vp=  435  m/s  measured  using  a neu- 
tron wavelength  of  A = 1.0  nm  (vn=  400  m/s).  The  contour  levels 
correspond  to  absolute  intensities  of  0.08,  0.03,  0.008  and  0.003 
cm  \ respectively. 


scatter  isotropically.  This  Doppler  shift  can  be  used  to 
directly  measure  the  velocity  of  the  particles,  and,  in  a 
test  case,  this  velocity  was  within  2 % of  the  average 
velocity  derived  from  the  pressure  trace  measurements 
[2]. 

Rather  than  simply  examine  nanodroplets  for  which 
segregation  of  the  components  is  severe,  we  use  different 
materials  to  observe  the  transition  between  well-mixed 
and  fully  segregated  droplets.  Thus,  our  experiments 
include  binary  mixtures  of  H,0,  D,0,  ethanol,  and  n- 
butanol,  or  its  fully  deuterated  analog  d-butanol.  Almost 


RESEARCH  HIGHLIGHTS 


B.  E.  Wyslouzil,  C.  H.  Heath, 

U.  M.  Dieregsweiler, 
and  K.  A.  Streletzky 

Department  of  Chemical  Engineering 
Worcester  Polytechnic  Institute 
Worcester,  MA  01609 


R.  Strey  and  J.  Wolk 

Institut  fur  Physikalische  Chemie 
Universitat  zu  Koln 
Cologne,  Germany 

G.  Wilemski 

Physics  Department 
University  of  Missouri  - Rolla 
Rolla,  MO  65409 


J.  Barker 

NIST  Center  for  Neutron  Research 
National  Institute  of  Standards  and  Technology 
Gaithersburg,  MD  20899-8562 


degenerate  mixtures,  such  as  D,0  - H,0,  appear  to  form 
well-mixed  droplets  [4].  In  contrast,  binary  nanodroplets 
containing  H,0  (or  D,0)  and  a low  molar  concentration 
of  d-butanol  (or  h-butanol)  should  consist  of  a water-rich 
core  and  an  alcohol-rich  shell.  Our  recent  SANS  experi- 
ments with  H,0  - d-butanol  nanodroplets  clearly  demon- 
strate for  the  first  time  the  existence  of  the  alcohol-rich 
shell.  As  illustrated  in  Fig.  3,  for  these  droplets  the 
scattered  intensity  falls  off  as  q'2  in  the  high  q region,  as 
is  characteristic  for  scattering  by  a shell.  Furthermore,  the 
feature  in  the  region  0.4  nnr'  < q < 0.5  nnr1  cannot  be 
reproduced  by  a well-mixed  droplet  model  that  matches 
the  low  q data.  Also  shown  in  Fig.  3 is  the  scattering 
spectrum  from  the  complementary  D,0  - h-butanol 
experiment.  As  expected,  the  scattered  intensity  is  much 
higher,  and  the  signal  intensity  decreases  as  q~4  in  the 
high  q region.  In  this  case,  the  D,0-rich  core  dominates 
the  scattering. 

In  summary,  SANS  provides  us  with  a powerful  new 
way  to  study  the  properties  of  nanometer  sized  liquid 
droplets  in  the  environment  in  which  they  form.  To  date, 
it  is  the  only  technique  to  directly  probe  the  microstruc- 
ture of  aerosol  nanodroplets.  Combined  with  pressure 
trace  measurements  and  modeling,  SANS  provides 
information  critical  to  our  understanding  of  droplet 
formation  and  growth  in  the  nanometer  size  regime. 
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Fig.  3.  Scattering  from  a D20  - h-butanol  aerosol  is  compared  to 
that  from  a H20  - d-butanol  aerosol  formed  under  identical  con- 
ditions in  the  nozzle.  Both  aerosols  contain  * 6 mol  % butanol. 
In  the  high  q region,  the  intensity  falls  off  as  qA  for  the  D20-rich 
droplets,  but  only  as  q2  for  the  H20-rich  droplets.  Our  prelimi- 
nary modeling  shows  that  the  feature  in  the  region  0.4  nm  1 < q < 
0.5  nm  1 of  the  H20  - d-butanol  spectrum  cannot  be  reproduced 
by  a well  mixed  droplet  model.  SDD  is  the  sample-to-detector 
distance. 
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Nanoporous  Ultra-Low  Dielectric  Constant 

Materials 


The  next  generation  of  interlayer  dielectric  materials 
for  microelectronics  must  have  an  ultra-low  dielectric 
constant  of  less  than  2.0  to  meet  the  National  Technology > 
Roadmap  for  Semiconductors  [ 1 ].  In  addition,  the 
material  must  have  excellent  thermal  stability  and 
mechanical  properties  [2.3],  One  route  is  to  introduce 
nanometer  sized  pores  into  a material  with  an  initially 
low  dielectric  constant  (k  ~ 2 to  3).  The  presence  of  the 
pores  lowers  the  dielectric  constant  of  the  nanoporous 
film  as  the  dielectric  constant  of  air  is  1.0.  One  class  of 
materials  that  is  receiving  significant  attention  is  based 
on  poly(methylsilsesquioxane)  (PMSSQ)  which  has  a 
dielectric  constant  between  2.7  and  2.8  (lower  than 
conventional  SiO.)  and  good  thermal  stability  up  to 
500  °C  [4,  5].  The  dielectric  constant  is  a function  of  the 
level  of  porosity  and  detailed  knowledge  of  the  morphol- 
ogy is  required  for  implementation  in  electronic  packaging. 
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porogen  material 
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precursor  material 
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Fig.  1 . Schematic  showing  the  generation  of  nanoporous  PMSSQ 
films. 

A schematic  diagram  of  the  synthesis  process  is 
shown  in  Fig.  1 . The  materials  are  based  on  an  inorganic 
matrix  of  PMSSQ  containing  a porogen  (i.e.  a labile  pore 
generating  material)  based  on  a poly( methylmethacrylate) 
copolymer  (coPMMA)  as  a template  for  creating  the 
pores  [6.7],  A film  is  prepared  by  dissolving  the  PMSSQ 
matrix  precursor  material  and  the  coPMMA  porogen  in  a 
common  solvent  and  spin  casting  onto  a silicon  substrate. 

The  film  is  then  slowly  ramped  up  in  temperature  to 
450  °C.  During  this  temperature  ramp,  the  matrix  film 


fully  cures  by  about  225  °C.  resulting  in  microphase 
separation  of  the  coPMMA  polymer.  Above  about  350  °C 
the  coPMMA  polymer  degrades  leaving  nanometer  sized 
pores.  Neutron  scattering  was  used  with  a deuterated 
coPMMA  polymer  to  understand  the  morphological 
development  of  the  films  during  the  cure.  The  use  of  the 
deuterated  porogen  provides  contrast  for  neutron  scatter- 
ing in  the  hybrid  (polymer  containing)  system  that  is  not 
available  with  x-ray  scattering. 

In  situ  studies  of  the  curing  process  by  small  angle 
neutron  scattering  (SANS)  were  performed.  To  provide 
sufficient  scattering  signal,  four  or  more  samples  on  Si 
wafers  were  stacked  to  increase  the  scattering  volume. 

Figure  2 shows  a set  of  SANS  data  from  a PMSSQ 
thin  film  containing  20  % porogen  at  various  cure 
temperatures.  The  scattered  intensity  varies  as  a function 
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Fig.  2.  SANS  data  taken  during  cure  of  a nanoporous  PMSSQ 
thin  film. 


of  curing  temperature  due  to  morphology  changes 
occurring  in  the  sample.  For  the  as-spun  materials  there 
is  a shoulder  at  high  q from  which  a radius  of  gyration 
(R  ) of  12  A was  obtained  by  Guiltier  analysis.  After  heat 
treatment  at  225  °C  the  shoulder  in  l(q)  becomes  more 
pronounced  and  shifts  to  about  q = 0.06  A'1.  The  shift  to 
smaller  q indicates  a coarsening  of  the  microphase 
separated  porogen  domain  structure.  A radius  of  gyration 
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Fig.  3.  Above,  neutron  reflectivity  data  (symbols)  and  best  fit 
(lines)  calculated  from  the  scattering  length  density  profiles 
shown  below  for  a nanoporous  PMSSQ  thin  film  under  various 
conditions. 

of  R = 23  A for  the  domains  was  obtained  for  the  cured 

g 

materials  . The  scattering  curves  do  not  change  signifi- 
cantly between  225  °C  and  300  °C,  indicating  that  the 
morphology  remains  fixed  until  the  degradation  of  the 
polymer  at  450  °C.  The  dramatic  decrease  in  intensity  of 


the  SANS  curve  at  450  °C  is  due  to  a loss  of  neutron 
scattering  length  density  contrast  because  of  the  degrada- 
tion of  the  porogen  and  the  formation  of  pores  in  the 
film. 

To  provide  neutron  contrast  between  the  pores  and  the 
matrix,  the  fully  cured  sample  was  exposed  to  deuterated 
toluene.  These  data  are  also  shown  in  Fig.  2.  The  scatter- 
ing is  largely  recovered  after  exposure  to  deuterated 
toluene,  indicating  that  the  porous  structure  is  maintained 
with  no  collapse  after  degradation  of  the  polymer. 

Figure  3 shows  a set  of  neutron  reflectivity  profiles 
and  scattering  length  density  (SLD)  profiles  of  a 
nanoporous  PMSSQ  thin  film  prepared  with  20  % 
porogen  and  deposited  on  a silicon  substrate.  The  film 
was  exposed  to  deuterated  toluene  liquid  and  vapor  to 
examine  solvent  swelling.  The  film  thickness  is  observed 
to  increase  upon  exposure  to  both  liquid  and  vapor.  A 
higher  porosity  is  observed  in  a 50  A region  next  to  the 
silicon  surface.  This  is  probably  due  to  porogen  surface 
segregation  at  the  substrate  interface  and  is  in  agreement 
with  results  from  transmission  electron  microscopy. 
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Neutron  Reflectivity  Studies  of  Surfactants  at 

Electrode  Surfaces 


Knowledge  of  the  forces  that  control  the  assembly  of 
surfactant  molecules  at  the  solid-liquid  interface  is  vital 
for  traditional  fields  such  as  detergency,  flotation,  oil 
recovery  and  tribology  1 1 ].  Thin  organic  films  deposited 
at  solid  surfaces  also  find  application  in  the  fabrication  of 
optoelectronic  devices,  sensors,  biosensors  and  chemi- 
cally modified  electrodes  [2],  It  has  long  been  established 
that  the  assembly  of  surfactants  at  the  solid-liquid 
interface  depends  on  the  charge  at  the  solid  surface  [3,4], 
For  example,  the  spreading  of  vesicles  into  a phospho- 
lipid bi  layer  requires  that  surface  of  the  solid  support  be 
negatively  charged  and  hydrophilic. 

However,  the  present  understanding  of  the  role  played 
by  charge  on  the  interaction  of  a surfactant  molecule  with 
the  electrified  solid  surface  is  far  from  being  complete. 
Electrochemistry  provides  a unique  opportunity  to  study 
the  effect  of  the  charge  on  the  properties  of  amphiphilic 
and  ionic  surfactants  at  the  charged  solid-liquid  interface 
[5],  When  an  organic  film  is  deposited  on  a gold  elec- 
trode the  charge  density  at  the  metal  surface  may  be 
varied,  from  about  30  pC/cnr  to  about  40  pC/cin2.  This 
magnitude  of  charge  generates  fields  on  the  order  of  10'° 
V/m.  Such  a field  interacts  with  polar  molecules  in  the 
membrane.  By  changing  the  sign  of  the  charge  one  can 
use  attractive  or  repulsive  forces.  In  this  manner,  by 
turning  a knob  on  a control  instrument  one  can  force 
phase  transitions  in  the  film  of  organic  molecules  or  force 
surfactants  to  desorb  or  re-adsorb  on  the  surface. 

We  have  recently  employed  electrochemical  tech- 
niques, atomic  force  microscopy  and  neutron  reflectivity 
to  study  the  field  driven  transformations  of  thin  films 
formed  by  a model  anionic  surfactant  sodium  dodecyl 
sulfate  (SDS)  at  the  surface  of  a gold  electrode  [6]. 

Figure  1 shows  how  the  surface  concentration  of  SDS  at 
the  gold  electrode  surface  changes  with  the  electrode 
potential.  A convenient  way  to  interpret  these  data  is  to 
treat  the  electrode  potential  as  an  operational  variable 
that  can  be  easily  adjusted  using  a control  instrument. 

Figure  1 shows  that  the  character  of  SDS  adsorption 
is  strongly  influenced  by  the  charge  on  the  metal.  At 
sufficiently  negative  potentials  SDS  molecules  are  totally 
desorbed  from  the  electrode  surface.  At  moderate 


A (mV  vs  sc E) 


Fig.  1.  Three  dimensional  plots  of  the  surface  concentration  of 
SDS  as  a function  of  electrode  potential  measured  versus  the 
calomel  reference  electrode  and  the  logarithm  of  the  bulk  SDS 
concentration. 

negative  charge  densities  SDS  forms  a film  characterized 
by  a limiting  surface  concentration  4.0  x I0i0  mol  cnr:. 
When  the  metal  surface  is  positively  charged  the  surface 
concentration  of  SDS  increases  to  8.1  x 10  i0  mol  cm  -. 

Neutron  reflectivity  was  employed  to  determine  the 
structure  of  the  film  formed  by  SDS  at  different  charge 
densities  at  the  gold  surface.  Neutron  reflectivity  experi- 
ments were  carried  out  on  the  NG-7  reflectometer.  Thin 
layers  of  chromium  («  20  A)  and  gold  (~  80  A)  were 
sputtered  onto  the  crystal  quartz  substrate.  After  cleaning, 
the  crystal  was  mounted  on  a specially  constructed  Teflon 
cell  1 7],  The  cell  had  ports  for  the  counter  (gold  foil)  and 
reference  electrodes  (Ag/AgCl,  E « - 40  mV  versus  SCE). 
D,0  (99.9  %)  was  used  as  a solvent  in  reflectivity 
studies. 

Figure  2a  shows  the  neutron  reflectivity  data  deter- 
mined for  SDS  adsorption  at  various  electrode  potentials, 
and  Fig.  2b  shows  the  scattering  length  density  profiles 
calculated  from  the  reflectivity  curves.  The  neutron 
reflectivity  data  are  consistent  with  electrochemical 
measurements.  They  show  that  at  very  negative  potentials 
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Fig.  2.  (a)  Normalized  neutron  reflectivity  curves  for  a Au/Cr- 
coated  quartz  substrate  in  16  mM  SDS  in  D20.  (b)  Scattering 
length  density  profiles  of  the  interface  as  determined  from  the 
fitting  procedure. 

the  gold  solution  interface  is  free  from  hydrogenated 
species.  When  the  potential  increases  the  film  of  hydro- 
genated species  appears  at  the  electrode  surface.  The 
thickness  of  this  film  increases  and  the  scattering  length 
density  progressively  decreases  with  increasing  potential. 
When  combined  with  the  results  of  electrochemical 
measurements  and  atomic  force  imaging,  the  neutron 
reflectivity  data  allow  the  determination  of  the  structure 
of  the  film  formed  at  different  charge  densities.  At  small 
or  moderate  negative  charge  densities,  SDS  molecules 
form  a hemimicellar  film  that  consists  of  hemicylindrical 
stripes,  as  first  observed  by  Mamie  [4].  The  packing  of 
SDS  molecules  in  a cross  section  of  that  hemicylinder  is 


shown  schematically  in  Fig.  3 (top).  At  positive  charge 
densities  the  hemimicellar  state  is  transformed  into  the 
interdigitated  bilayer  schematically  shown  in  Fig.  3 
(bottom). 

The  results  of  this  study  demonstrate  the  need  for  the 
use  of  neutron  reflectometry  to  explore  adsorption  and 
phase  transitions  in  films  of  surfactants  adsorbed  at  the 
solid-solution  interface.  Specifically,  they  show  that 
when  neutron  reflectivity  measurements  are  combined 
with  electrochemical  studies  and  atomic  force  micros- 
copy, they  provide  unique  opportunities  to  explore 
different  stages  involved  in  the  interaction  of  surfactants 
with  solid  surfaces. 
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Fig.  3.  Models  of  SDS  adsorption  at  the  Au-  solution  interface:  a 
cross  section  of  hemicylindrical  aggregates  observed  at  moder- 
ately negative  charge  densities  (top);  interdigitated  bilayer 
observed  at  positive  charge  densities  (bottom). 
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Magnetic  Semiconductor  Superlattices 


Currently  a great  deal  of  attention  is  being  focused  on 
spintronics,  a new  area  of  solid-state  electronics.  In 
spintronics  not  only  the  current  magnitude  but  also  its 
spin  state  is  controlled.  Spin  valves  and  spin  injectors  are 
the  first  practical  applications  of  spintronics.  Further 
progress  in  developing  new  devices  hinges  critically  on 
the  availability  of  suitable  materials.  Such  materials  need 
to  be  “good"  semiconductors,  easy  to  integrate  in  typical 
integrated  circuits,  and  their  electronic  properties  should 
exhibit  strong  sensitivity  to  the  carrier’s  spin,  ferromag- 
netism being  an  especially  desirable  property. 

EuS  is  one  of  the  very  few  natural  ferromagnetic  (FM) 
semiconductors.  Since  it  becomes  FM  at  a low  tempera- 
ture (T  - 16.6  K)  it  is  an  unlikely  choice  for  applica- 
tions. However,  studying  the  properties  of 
heterostructures  made  on  its  base  may  give  an  important 
insight  into  fundamental  processes  taking  place  in  all 
classes  of  materials  under  consideration. 

GaMnAs  is  a man-made  FM  semiconductor.  It  is  an 
example  of  a diluted  magnetic  semiconductor  (DMS)  in 
which  a fraction  of  nonmagnetic  cations  (Ga)  is  substi- 
tuted with  magnetic  ions  (Mn).  Such  a material  can 
readily  be  incorporated  into  modern  GaAs-based  semi- 
conductor devices.  Its  Tc  is  still  below  room  temperature, 
but  this  limitation  may  be  lifted  in  other  materials  of  this 
class  [ I ]. 

Interlayer  exchange  coupling  (IEC)  in  superlattices 
(SL),  composed  of  ferromagnetic  and  nonmagnetic 
layers,  is  a crucial  element  of  all  spin-valve  type  devices 
that  utilize  the  giant  magnetoresistance  effect.  In  metallic 
SLs  currently  being  used,  conduction  electrons  transfer 
the  interlayer  interactions  through  nonmagnetic  spacers 
[2].  Here  we  address  the  question  whether  IEC  phenom- 
ena are  possible  in  all-semiconductor  superlattices,  like 
EuS/PbS  and  GaMnAs/GaAs,  where  the  carrier  concen- 
trations are  many  orders  of  magnitude  lower  than  in 
metals. 

The  nonmagnetic  spacer  in  EuS/PbS  SLs  is  a narrow 
gap  semiconductor  with  electron  concentration  of  the 
order  of  1017  cm  3 to  1018  cm  3.  For  thin  PbS  layers  (dpbS< 
70  A)  neutron  reflectivity  spectra,  shown  in  Fig.  1,  have 
revealed  a pronounced  maximum  of  magnetic  origin  at 


the  position  corresponding  to  the  doubled  structural  SL 
periodicity,  thus  indicating  the  existence  of  antiferromag- 
netic (AFM)  interlayer  arrangements  [3], 

For  much  thicker  PbS  spacers  (dpbS>  120  A)  the  only 
magnetic  peaks  visible  in  the  reflectivity  profiles,  see  e.g. 
Fig.  2,  coincide  with  the  chemical  ones,  thus  leading  to 
the  conclusion  that  the  magnetization  vectors  in  adjacent 
EuS  layers  are  parallel  which  indicates  FM  IEC. 
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Fig.  1.  Unpoiarized  neutron  reflectivity  spectra  for  EuS/PbS  SL 
with  thin  (23  A)  PbS  spacer.  Antiferromagnetic  interlayer  ex- 
change coupling  below  Tc  and  at  zero  external  field  is  clearly 
visible  (blue  curve).  Applying  a strong  enough  magnetic  field 
(185  G in  this  case)  parallel  to  the  SL  surface  forces  all  the  EuS 
layer’s  magnetizations  to  ferromagnetic  configuration  (green 
curve).  Above  Tc  the  system  is  nonmagnetic,  the  only  Bragg  peak 
comes  from  the  chemical  SL  periodicity. 

In  the  intermediate  PbS  thickness  range  (70  A < dphS< 
120  A),  both  types,  AFM  as  well  as  FM,  peaks  are 
present.  Polarized  neutron  analysis  of  these  maxima  gives 
evidence  that  the  magnetization  vectors  of  adjacent  EuS 
layers  are  not  collinear.  Hence,  the  IEC  found  in  EuS/PbS 
SLs  has  an  oscillatory  character  similar  to  that  occurring 
in  metallic  SLs,  although  the  oscillation  period  is  much 
longer  than  the  one  in  metallic  systems. 

In  order  to  confirm  that  the  free  carriers,  present  in 
the  PbS  layer  in  such  a small  amount,  are  the  cause  of  the 
observed  oscillatory  IEC,  a series  of  analogous  measure  - 
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Fig.  2.  The  sample  with  thick  (135  A)  PbS  layers  is  almost 
ferromagnetically  coupled.  Application  of  an  external  magnetic 
field  enhances  the  FM  Bragg  peaks  and  lowers  the  intensity  be- 
tween them  (in  AFM  peak  position). 


ments  have  been  carried  out  on  EuS/YbSe  SLs.  The 
structure  and  lattice  constant  of  YbSe  are  the  same  as 
those  of  PbS.  In  contrast  to  PbS,  YbSe  is  a semi-insulator 
with  a negligible  carrier  concentration.  Neutron 
reflectivity  profiles  have  shown  no  evidence  of  any 
interlayer  coupling  in  the  all  investigated  samples.  That 
finding,  together  with  the  oscillatory  character  of  cou- 
pling in  SLs  with  PbS  spacer,  strongly  points  to  the 
leading  role  of  PbS  free  electrons  in  providing  the 
necessary  IEC  mechanism,  similar  to  that  discovered  in 
metallic  multilayers. 

Ferromagnetic  ordering  in  GaMnAs  is  carrier  (holes) 
induced;  the  Mn  atoms,  apart  from  being  the  magnetic 
element  in  the  system,  act  also  as  acceptors  providing  the 
holes  responsible  for  transferring  exchange  interactions 
between  them.  The  details  of  the  magnetic  ordering,  in 
particular  its  range,  are  still  being  disputed. 

To  address  the  latter  issue,  polarized  neutron  reflecto- 
metry  has  been  performed  on  a number  of  GaMnAs/ 
GaAs  superlattices.  Figure  3 shows  an  example  of  the 
reflectivity  profile  obtained  in  the  vicinity  of  the  first  SL 
Bragg  peak  for  one  of  the  samples.  The  very  presence  of 


the  magnetic  contribution  to  the  structural  SL  Bragg  peak 
is  a strong  confirmation  of  the  FM  IEC  between  consecu- 
tive GaMnAs  layers.  The  absence  of  any  spin-flip 
scattering  shows  that  the  sample  is  in  a one-domain  state, 
i.e.,  the  FM  ordering  in  GaMnAs  is  long  range,  and  the 
sample  is  spontaneously  saturated.  The  peak  in  (— ) cross 
section,  and  its  absence  in  the  (++),  is  proof  that  the 
magnetization  is  directed  oppositely  to  the  external 
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Fig.  3.  Polarized  neutron  reflectivity  spectra  for  GaMnAs/GaAs 
superlattice. 


magnetic  guide  field,  hence  the  long  range  ordering  has 
formed  spontaneously,  without  the  influence  of  the 
external  field.  More  details  can  be  found  in  [4], 
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Pinpointing  Chiral  Structures  with  Front/Back 

Polarized  Neutron  Reflectometry 


We  have  developed  a new  method  of  using  polarized 
neutron  reflectometry  (PNR)  to  extract  the  structure  of 
buried  magnetic  spirals  in  magnetic  films.  This  technique 
improves  upon  earlier  methods  by  being  particularly 
sensitive  to  the  presence  of  magnetic  twists  vis-a-vis 
structures  in  which  the  magnetization  direction  does  not 
vary  appreciably.  Tracking  the  formation  and  growth  of 
twists  may  solve  a number  of  puzzles  that  hamper  the 
development  of  magnetic  thin  film  devices. 

In  collaboration  with  IBM  scientists,  we  have  applied 
the  technique  to  a thin-film  exchange-spring  magnet  and 
confirmed  that  the  results  may  violate  the  current  theory 
regarding  the  behavior  of  such  magnets.  It  has  been 
predicted  that  exchange-spring  magnets,  comprised  of 
soft  and  hard  ferromagnets  in  close  proximity,  are  a 
composite  that  has  a strong  moment  and  does  not  readily 
demagnetize  [1 ).  Therefore,  exchange-spring  magnets 
should  give  industry  the  ability  to  make  much  smaller 
permanent  magnets  for  use  in  the  magnetic  recording 
devices,  and  elsewhere.  As  a side  effect,  when  a small 
external  magnetic  field  is  opposed  to  that  of  the  magnet, 
the  portion  of  the  soft  ferromagnet  farthest  from  the  hard 
ferromagnet  may  twist  into  alignment  with  the  field. 

When  the  field  is  removed,  the  soft  ferromagnet  untwists. 
The  film  provided  by  IBM  consists  of  the  hard  ferromag- 
net Fe55Pt4.  topped  by  the  soft  ferromagnet  Ni80Fe,0  [2], 

Figure  1 shows  a simplified  diagram  of  the  behavior 
predicted  by  current  theories  [1  ].  A magnetic  field  of 
0.890  T,  provided  by  an  electromagnet,  is  sufficient  to 
align  both  the  soft  and  the  hard  layers  of  our  exchange- 
spring magnet,  as  shown  on  the  left.  When  a modest 
reverse  field  (on  the  order  of  0.025  T)  is  applied  to  the 
exchange  spring  magnet,  only  the  top  of  the  soft  layer 
will  realign  with  the  magnetic  field.  The  hard  layer 
remains  pinned  in  the  original  direction,  and  a continuous 
twist  is  induced  in  the  soft  layer,  as  the  direction  of 
magnetization  changes  smoothly  between  the  reverse 
field  direction  to  the  aligning  field  direction. 

Although  there  are  many  alternatives  to  PNR  to 
measure  the  magnetization,  typically  they  measure  only 
the  average  orientation  of  the  magnetic  spins,  and  cannot 
readily  distinguish  a spiral  from  a structure  in  which  all 


Hard  ferromagnet 


Fig.  1.  Model  for  field  behavior  of  exchange-spring  magnets.  On 
the  left  the  magnet  has  been  aligned  by  a large  external  mag- 
netic field.  On  the  right  a smaller  field  opposed  to  the  first  field 
causes  a twist  to  form  in  the  soft  ferromagnet,  while  the  hard 
ferromagnet  remains  aligned. 

the  spins  are  canted  with  respect  to  an  external  field. 

PNR  can  extract  the  depth-dependence  of  the  magnetic 
and  chemical  structure.  We  have  studied  the  sample  over 
a wide  range  of  external  magnetic  fields,  and  can  track 
the  development  of  the  spiral  with  field  [3], 

A PNR  experiment  begins  with  neutrons  whose 
magnetic  moments  are  aligned  parallel  (+)  or  opposite  (-) 
to  the  external  magnetic  field.  When  the  magnetization  of 
the  sample  is  perpendicular  to  this  magnetic  field,  the 
neutron  moment  precesses  as  it  interacts  with  the  sample. 
When  this  happens  the  spin-flip  (SF)  reflectivities  R + ~ 
and  R~ + are  strong.  If  the  magnetization  of  the  sample  is 
parallel  to  the  external  magnetic  field,  no  precession 
occurs,  but  the  non-spin-flip  (NSF)  reflectivities  R + + and 
R will  differ.  The  NSF  reflectivities  also  provide 
information  about  the  chemical  structure  of  the  film. 

Our  new  modification  of  the  PNR  method  greatly 
enhances  the  contrast  between  col  I inear  and  certain 
noncollinear  magnetic  structures  [4],  We  first  measure 
the  reflectivity  with  neutrons  glancing  off  the  front 
surface  of  the  material,  and  then  repeat  with  neutrons 
glancing  off  the  back  surface.  The  experiment  is  akin  to 
holding  the  plane  of  the  film  up  to  a “magnetic  mirror”  to 
see  whether  the  mirror  image  is  the  same  as  the  original 
structure.  In  a collinear  structure,  all  the  spins  are  aligned 
along  a common  direction,  and  the  mirror  image  is  very 
much  like  the  original  structure.  But  the  mirror  image  of 
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Fig.  2.  Reflectivity  of  a Ni80Fe20|Fe55Pt45  bilayer.  The  front  plotted  against  the  right  ordinate  axis.  The  NSF  reflectivities 

reflectivity  is  plotted  on  the  right  while  the  back  reflectivity  R.t  gnd  R are  plotted  against  the  ,eft  ordinate  axis, 

is  plotted  on  the  left.  The  SF  reflectivities  R*  and  R * are 


a magnetic  twist  to  the  right  is  a magnetic  twist  to  the 
left.  Therefore,  if  the  front  and  back  reflectivities  are 
significantly  different,  we  can  deduce  the  presence  of  a 
spiral.  Fitting  the  data  confirms  the  spiral's  existence. 

Figure  2 shows  data  collected  at  0.026  T after  aligning 
in  -0.89  T.  Fits  to  the  data  are  shown  as  solid  lines.  The 
data  from  the  front  reflectivity  are  shown  on  the  right, 
and  the  data  from  the  back  reflectivity  are  shown  on  the 
left.  The  spin-flip  (SF)  reflectivities  R ¥~  and  R plotted 
against  the  right-hand  axis,  have  been  shifted  relative  to 
the  NSF  reflectivities  R++  and  R that  are  plotted  against 
the  left  axis.  At  q = 0.2  nnr1,  there  is  a splitting  in  the 
front  NSF  reflectivity  that  is  much  more  pronounced  than 
that  of  the  back  reflectivity  at  the  same  q.  This  is  a 
hallmark  of  the  spiral  structure. 

Figure  3 shows  the  magnetic  structure  that  gives  the 
excellent  fit  to  the  data  plotted  in  Fig.  2.  The  location  of 
the  hard/soft  interface  is  marked  in  Fig.  3.  Surprisingly, 
we  discover  the  spiral  invades  the  hard  ferromagnet  even 
at  extremely  low  fields.  Current  theory  predicts  that  when 
this  occurs,  the  soft  ferromagnet  will  not  be  able  to 
untwist  fully.  Yet,  other  magnetic  studies  show  that  our 
exchange-spring  magnet  does  untwist  when  this  field  is 
removed.  Thus,  our  PNR  measurements  have  identified  a 
shortcoming  of  current  theory. 

With  this  new  technique,  NIST  is  now  able  to  better 
characterize  the  magnetic  properties  of  thin  films,  which 


can  improve  the  capability  and  reliability  of  industrial 
devices  for  magnetic  recording  and  sensing. 


Fig.  3.  Fitted  magnetization  of  the  data  presented  in  Fig.  2.  The 
front  of  the  sample  is  at  a depth  of  0 nm  and  the  back  is  at  a 
depth  of  70  nm.  The  red  curve  is  a projection  of  the  magnetic 
structure  into  the  plane  of  the  front  surface. 
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SF  Refectivity 


The  Nature  of  Vibrational  Softening  in 


a-Uranium 

The  standard  textbook  explanation  for  phonon 
softening  with  increasing  temperature  in  a single  phase  is 
that  the  potentials  are  not  perfectly  harmonic,  but  it  has 
been  suggested  that  phonon  softening  can  also  occur  if 
the  potential  itself  can  change  with  temperature,  while 
remaining  harmonic.  For  example,  a large  softening  of 

uranium  at  the  four  highest  temperatures  scaled  so  that 
the  lowest  point  is  at  the  harmonic  energy  kHT/2.  The 
effects  of  anharmonicity  on  <U>  would  be  evident  as 
nonlinearity  in  the  plot  of  <U>  vs.  temperature.  For 
comparison,  attempts  were  made  to  calculate  the 

Ihe  vibrational  properties  of  a-uranium  has  been  ob- 
served 1 1]  that  cannot  be  explained  in  terms  of  the 
anharmonic  lattice  contribution  alone.  In  the  present 
work  a harmonic  contribution  to  the  phonon  softening  is 
made  evident  by  treating  inelastic  neutron  scattering 
spectra  as  an  expansion  of  the  vibrational  power  spec- 
trum of  the  atomic  motion. 

Earlier  diffractometry  measurements  by  Lawson  et  al. 
[2]  suggested  that  the  Debye  temperature  decreased  by 
40  % in  going  from  300  K to  the  (3-phase  transition 
temperature  at  940  K.  In  terms  of  entropy  this  corre- 
sponds to  an  additional  AS  = -3kjn(0.6)  = 1.5  kB/atom. 
The  usual  thermodynamic  argument  is  that  this  increase 
in  vibrational  entropy  compensates  for  the  elastic  energy 
generated  by  thermal  expansion.  However,  from  heat 
capacity  data,  the  entropy  needed  to  compensate  for  the 
elastic  energy  is  nearly  an  order  of  magnitude  too  small 
at  0.16  kjatom. 

D 

All  experiments  used  uranium  powder.  High  tempera- 
ture measurements  were  made  using  the  Fermi-Chopper 
Spectrometer  (FCS)  at  the  NCNR.  Low  temperature 
measurements  were  performed  with  the  Low  Resolution 
Medium  Energy  Chopper  Spectrometer  (LRMECS)  at 
Argonne  National  Laboratory.  Figure  1 shows  the  phonon 
density  of  states  obtained  from  the  measured  spectra, 
corrected  for  the  incoherent  multiphonon  scattering  using 
a procedure  described  elsewhere  [3],  There  is  a redistri- 
bution of  intensity  in  the  main  features  at  « 8 meV  and 
« 12  meV,  with  the  higher  energy  peak  gaining  extra 
weight  with  increasing  temperature.  These  features  also 
show  an  overall  softening  of  around  one  meV  for  every 
200  K increase  in  temperature. 

The  Q-summed  one  phonon  scattering  function  was 
used  to  calculate  a quantity  proportional  to  the  square  of 
the  power  spectrum  and  hence  to  the  average  potential 
energy  per  oscillator,  <U>.  In  Fig.  2 we  show  <U>  for  a- 


5 10  15 

Energy  Transfer  (meV) 

Fig.  1 . The  phonon  density  of  states  of  uranium.  Data  from  300  K 
and  above  were  obtained  from  spectra  acquired  with  FCS.  Data 
from  300  K and  below  were  measured  on  LRMECS. 
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temperature  dependence  of  the  potential  energy  of  Morse 
and  Lennard-Jones  potentials  with  the  appropriate 
vibrational  softening  for  uranium  [I],  The  result  shown 
in  Fig.  2 indicates  that,  in  the  high  temperature  limit,  the 
potential  energy  has  nearly  a linear  dependence  on 
temperature,  i.e.  the  phonon  softening  in  a-uranium 
occurs  while  the  potential  remains  primarily  harmonic. 
Evidently  the  interatomic  force  constants  are  decreasing 


Fig.  2.  Vibrational  potential  energy  of  a-uranium  (o).  The  Lennard- 
Jones,  Morse_1  and  Morse_2  curves  were  calculated  from 
potentials  described  in  the  text.  The  Harmonic  curve  is  the  re- 
sult for  a harmonic  potential  in  the  classical  limit. 

with  increasing  temperature.  Since  the  force  constants 
originate  from  the  change  of  the  electronic  energy  with 
atom  displacements,  it  must  be  that  thermal  excitations 
of  the  electronic  states  are  altering  the  force  constants. 

The  phonon  density  of  states  of  the  three  solid  state 
phases  of  uranium,  orthorhombic  (a),  tetragonal  ( (3 ) and 
body  centered  cubic  (y)  are  compared  in  the  top  panel  of 
Fig.  1 . The  y-uranium  phonon  density  of  states  was 
statistically  identical  at  1 1 13  K and  1213  K.  Evidently, 
the  thermal  softening  mechanism  seen  in  a-phase  does 
not  operate  in  the  y-phase.  The  (3-phase  was  not  stable 


over  a wide  enough  temperature  range  to  obtain  a reliable 
temperature  dependence.  The  phonon  softening  between 
each  phase  accounted  for  vibrational  entropy  changes  of 
0.15  =t  0.01  Ay/atom  and  0.20  ± 0.0 1 Adatom  for  the  a^(3 
and  P— »y  transitions,  respectively.  Both  of  these  values 
make  up  only  about  35  % to  40  % of  the  total  entropy 
changes  predicted  from  latent  heat  measurements; 

(SP-Sa)lol  = 0.37  kj atom  and  0.55  kJdAom. 

The  remaining  60  % of  the  entropy  increases  must  be 
electronic  in  origin.  So  not  only  does  the  phonon  soften- 
ing disappear  in  the  high  temperature  y-phase,  but  also  it 
does  so  with  a large  increase  in  electronic  entropy. 

Electronic  band  structure  calculations  used  to  predict 
phonon  frequencies  are  based  on  the  assumption  that 
thermal  effects  can  be  neglected  when  compared  to 
volume  effects.  The  actinides,  however,  show  the  need 
for  more  sophisticated  treatments  of  the  role  of  tempera- 
ture on  interatomic  interactions. 
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Dynamics  of  Native  and  Partially  Unfolded  Proteins 
by  inelastic  Neutron  Scattering  and  Molecular 

Dynamics  Simulations 


A molten  globule  (MG)  is  a partially  folded  protein 
that  is  compact  with  native-like  secondary  structure  and 
overall  backbone-folding  topology  but  which  lacks  the 
extensive  specific  side-chain  packing  interactions  of  the 
native-like  structure.  Such  partially  folded  states  have 
been  shown  to  resemble  the  kinetic  folding  intermediate 
states  along  the  protein-folding  pathway.  In  addition, 
partially  folded  proteins  are  important  in  cell  functioning. 
A full  understanding  of  the  mechanism  of  protein  folding 
requires  the  knowledge  of  the  structures,  relative  energet- 
ics, and  dynamics  of  the  species  populating  the  folding 
pathway.  Incoherent  inelastic  neutron  scattering  (INS), 
which  makes  use  of  the  large  incoherent  cross  section  of 
hydrogen  nuclei,  is  a technique  well  suited  to  the  study  of 
protein  internal  molecular  motion  on  the  picosecond  time 
scale.  Molecular  dynamics  (MD)  simulations  on  the  other 
hand  are  a potentially  valuable  tool  for  interpreting 
neutron  data  on  proteins  [1,2].  and  may  provide  atomic 
level  description  of  the  motions  taking  place. 

In  this  study,  we  have  explored  the  dynamics  of 
bovine  a-lactalbumin,  a calcium-binding  protein,  in  the 
native  and  molten  globule  states,  using  both  techniques. 
INS  data  were  collected  on  the  Disk-Chopper  Spectrom- 
eter at  a 32  peV  resolution  for  native  and  molten  globu- 
lar bovine  a-lactalbumin.  The  simulations  results  were 
generated  from  a 1 ns  constant  temperature  and  pressure 
MD  trajectories.  The  model  systems  consisted  each  of  a 
single  protein  monomer,  immersed  in  a large  water  box 
(60  A x 50  A x 80  A).  For  the  native  state,  the  a-lactalbu- 
min  configuration  was  initiated  from  the  crystal  structure. 
The  model  molten  globule  state  was  generated  by 
partially  unfolding  the  protein  at  500  K at  atmospheric 
pressure  until  the  radius  of  gyration  reached  the  experi- 
mental value  determined  for  the  molten  globule  state 
( 10  % expansion).  A 1 ns  trajectory  was  then  generated  at 
300  K for  data  analysis  and  comparison  with  experiment. 

Figure  1 shows  the  structure  of  the  a-lactalbumin  in 
the  native  and  molten  globule  states  obtained  from  the 
MD  simulations  in  solution  at  300  K.  The  overall  struc- 
tural change,  manifested  by  the  conservation  of  the 
helical  a-domain  and  the  loss  of  the  (3  sheets  in  the  P- 
domain,  is  consistent  with  NMR  data  [3],  We  therefore 


regard  this  structure  as  a reasonable  model  of  a member 
of  the  MG  ensemble  for  comparison  with  the  INS  data. 

In  Fig.  2 we  report  the  incoherent  structure  factors 
S(Q,co)  as  a functions  of  energy  transfer  Hco  for  selected 
momentum  transfer  Q , measured  and  calculated  from  the 
MD  simulations  trajectories.  In  both  systems  the  agree- 
ment with  experiment  is  remarkable,  assessing  the  ability 
of  the  simulations  to  reproduce  the  picosecond  dynamics 
of  the  protein.  In  agreement  with  previous  INS  data 
collected  at  lower  resolution  [4],  the  molten  globule  has  a 
broader  quasi-elastic  spectra  compared  to  the  native  state. 


Fig.  1.  a-lactalbumin  structure  in  the  native  (top)  and  molten  glob- 
ule (bottom)  state  from  MD  simulations.  Note  the  loss  of  the 
p-domain  structure,  with  most  the  a domain  (helix)  being  con- 
served. 

which  indicates  that  the  protein  atoms  are  more  mobile  in 
the  MG  state.  The  additional  motion  is  large  enough  to  be 
detected  at  the  length  scale  probed  by  the  DCS  spectrom- 
eter, /. e. , in  the  range  3 A to  60  A. 

To  investigate  the  relationship  between  the  structure 
of  the  protein  and  the  motion  of  its  side  chains,  we  have 
used  the  MD  trajectories  to  follow  the  motion  of  these 
side  chains  quantitatively.  Figure  3 reports  the  fluctua- 
tions per  residue  along  the  backbone  for  the  native  and 
the  molten  globule  states,  as  a function  of  residue.  These 
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Fig.  2.  Incoherent  structure  factor  S(Q,m)  at  selected  wavevectors 
Q,  measured  and  calculated  from  the  simulation  for  the  native 
(top)  and  molten  globule  (bottom)  states. 

correspond  to  the  100  ps  root-mean-  square  amplitudes  of 
motion  of  the  backbone  atoms,  averaged  over  the  simula- 
tion trajectories.  First,  for  the  native  state,  the  backbone 
fluctuations  appear  to  correlate  with  the  secondary 
structure  of  the  protein,  showing  much  smaller  ampli- 
tudes in  the  a-helix  and  [3-sheet  regions.  For  the  molten 
globule  state,  a significant  increase  in  the  protein  fluctua- 
tion amplitudes  is  observed,  in  agreement  with  the 
broadening  of  the  INS  spectra.  Figure  3 shows  clearly 
that  most  of  the  additional  motion  occurs  in  the  region  of 
the  protein  that  “unfolds”,  i.e.,  the  (3-domain.  The  results 
indicate  also  that  the  motion  of  the  side  chains  is  by  no 
means  homogeneous.  Differences  in  amplitudes  up  to  10 
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Fig.  3.  100  ps  mean  squared  fluctuations  of  the  backbone  atoms 
computed  from  the  MD  simulations  of  the  native  state  (black) 
and  the  molten  globule  state  (red). 

folds  are  observed  between  parts  of  the  proteins,  regard- 
less of  exposure  to  the  solvent. 

The  MD  study  of  the  dynamics  of  partially  folded 
states  of  the  protein  should  in  principle  be  extended  to 
explore  a large  portion  of  the  conformational  space,  i.e., 
be  extended  to  a large  ensemble  of  non-native  conform- 
ers.  At  any  rate,  the  present  results,  though  limited, 
combined  with  neutron  experiments  are  believed  to 
capture  the  essence  of  the  structural  and  dynamical 
changes  taking  place  as  the  protein  partially  unfolds. 

They  demonstrate  the  utility  of  the  MD  simulations  for 
qualitatively  elucidating  the  dynamical  behavior  of  native 
and  non-native  proteins  at  the  atomic  level. 
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NIST  CENTER  FOR  NEUTRON  RESEARCH 


Dynamics  in  Confined  Polymer  Films 


Material  properties  at  surfaces  and  interfaces  often 
deviate  from  their  bulk  values.  This  has  ramifications  in  a 
range  of  polymer  thin  film  applications  such  as  chemi- 
cally amplified  photoresists,  lubricants,  adhesives,  etc.  Of 
particular  interest  are  changes  in  the  glass  transition 
temperature,  T , and  the  viscoelastic  properties  that  are 
related  to  7\  A viscosity  increase  of  liquid  polymeric 
films,  like  those  used  for  hard  disc  lubricants,  could  have 
deleterious  consequences.  Likewise,  a T decrease  in  a 
thin  photoresist  film  could  enhance  photoacid  diffusion 
and  create  image  blur.  Clearly,  it  is  of  great  interest  to 
characterize  the  molecular  level  mobility  in  thin  polymer 
films. 

It  is  well  known  that  thermal  motions  lead  to  a 
decrease  in  the  intensity  of  the  elastically  scattered 
neutrons.  In  a harmonic  solid  this  is  related  to  the 
amplitude  of  the  motion  through  the  Debye- Waller  factor: 
I (Q)  ~ exp (-'/}Q2<ir>)  where  O is  the  scattering  vector 
and  < ir>  is  the  mean-square  displacement.  In  a logarith- 
mic plot  of  Imc  as  a function  of  Q2,  <ir>  is  proportional 
to  the  slope  of  a linear  fit  to  the  data.  This  approximation 
has  proven  very  useful  in  studying  dynamics  of  polymers, 
both  synthetic  [1 1 and  biological  [2],  despite  the  fact  that 
their  motions  are  typically  anharmonic. 

The  elastic  incoherent  neutron  scattering  intensities 
can  be  measured  as  a function  of  temperature  on  the 
NCNR's  High  Flux  Backscattering  Spectrometer  (HFBS), 
using  the  “Fixed  Window”  mode.  Polycarbonate  (PC) 
films,  seen  in  Fig.  1,  are  spin  cast  onto  Si  wafers,  sec- 
tioned into  small  rectangles,  and  stacked  into  an  A I 
sample  cell.  The  scattering  from  the  Si  substrates  and  Al 
cell  are  negligible  because  their  Bragg  peaks  are  beyond 
the  O range  of  the  instrument  (0.25  < Q (A1)  < 1 .75);  the 
scattering  is  dominated  by  the  large  incoherent  cross- 
section  of  the  amorphous  PC. 

Ideally  there  should  be  200  mg  of  PC  per  cell  to 
maximize  the  signal  and  minimize  multiple  scattering 
effects  (90  % transmission).  For  the  Si-supported  films 
we  are  limited  by  ( 1 0 to  0.5)  mg  of  PC  per  cell  and  thus  a 
very  weak  signal  (99  % transmission).  However,  the 
HFBS  is  sufficiently  sensitive,  as  seen  in  Fig.  2 where  the 
temperature  dependence  of  the  elastic  scattering  is 


HFBS  cell  w/  1014  A 
PC  film  on  Si  wafers 


Fig.  1.  A stack  of  1014  A PC  films  (notice  the  blue  tint)  on  thin 
slabs  of  Si  ready  to  go  into  the  HFBS  spectrometer. 

shown.  The  relative  magnitude  of  the  intensity  change  is 
diminished  with  decreasing  film  thickness,  a clear 
indication  of  reduced  mobility  in  the  thin  PC  films.  Given 
the  HFBS  energy  resolution  of  0.85  peV,  these  changes 
reflect  atomic  motions  at  200  MHz  or  faster. 
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Fig.  2.  The  normalized  elastic  intensities  in  the  HFBS  detector  at 
Q = 0.36  A1.  The  decrease  in  intensity  at  high  temperatures 
reflects  dynamics  in  the  thin  polymer  films. 

This  reduced  mobility  is  qualitatively  reflected  in  Fig. 

3 as  a strong  suppression  of  <u2>.  We  emphasize  qualita- 
tively because  <u2>  is  derived  from  the  harmonic 
approximation.  Nevertheless,  there  is  a clear  reduction  in 
the  amplitude  of  the  thermal  motions  both  above  and 
below  the  calorimetric  T as  the  film  thickness  ap- 
proaches the  rnis  end-to-end  length  of  PC  (R  = 160  A). 
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At  temperatures  below  200  K,  <u2>  is  small  and  the 
harmonic  approximation  is  reasonable.  In  this  regime  the 
thermal  evolution  of  <ir>  (i.e.,  slope)  defines  an  elastic 
force  constant  k.  The  inset  of  Fig.  3 displays  a dramatic 
stiffening  of  a:  with  confinement.  This  implies  that  the 
atoms  or  molecules  in  PC  become  strongly  caged  by  their 
environment  in  the  exceedingly  thin  films. 

Similar  <u2>  measurements  on  thin  films  of  poly- 
methyl methacrylate  (PMMA;  R = 915  A)  and  polyvinyl 
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Fig.  3.  <u2>  as  a function  of  temperature  for  the  PC  thin  films. 
The  inset  reveals  a stiffening  of  the  harmonic  (low  temperature) 
force  constant  k with  decreasing  film  thickness  h. 


chloride  (PVC;  R = 437  A)  are  shown  in  Figs.  4(a)  and 
4(b).  There  is  a suppression  of  <ir>  with  decreasing  film 
thickness  above  the  bulk  T that  is  qualitatively  consistent 
with  the  PC  films.  However,  below  T bulk  the  suppres- 
sion is  less  pronounced  in  PMMA,  and  not  observed  in 
the  PVC  films.  Careful  inspection  of  Figs.  3 and  4 
suggests  a relation  between  the  magnitude  of  <ir>  at  T 
in  the  bulk  and  the  degree  of  <u2>  suppression  in  the 
thin  glassy  films. 

Bulk  PC  has  the  largest  <ir>  at  T , indicating  exten- 
sive mobility.  In  fact,  PC  is  renowned  for  is  segmental 
mobility  in  the  glassy  state,  which  is  in  pail  related  to  its 
superior  impact  resistance  (bullet  proof  glass);  molecular 
mobility  helps  dissipate  strain  energy.  This  extensive 
mobility  is  also  susceptible  to  thin  film  confinement,  as  is 
also  seen  through  the  stiffening  of  the  elastic  constants 
from  1.9  N/m  in  bulk  PC  to  5.8  N/m  in  the  75  A film 
(stronger  molecular  caging).  At  the  other  extreme,  bulk 
PVC  has  the  smallest  <ir>  at  T , indicating  a lack  of 
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Fig.  4.  <u2>  as  a function  of  thickness  and  temperature  for  simi- 
lar PMMA  (a)  and  PVC  (b)  films. 


motions  that  can  be  affected  by  confinement.  This  is  also 
consistent  with  the  strong  intermolecular  caging  at  low 
temperatures;  both  bulk  PVC  and  the  220  A film  have 
large  elastic  force  constants  of  4.6  N/m.  This  lack  of 
mobility  reflects  the  observation  that  PVC  is  a brittle 
polymer  that  needs  high  levels  of  plasticization  for 
practical  use.  Only  when  PVC  is  taken  above  the  bulk  7), 
where  long-range  motions  dominate,  do  the  effects  of 
confinement  become  evident. 

In  short,  we  have  seen  two  different  responses  to  thin 
film  confinement.  While  there  is  a general  suppression  of 
<u2>  and  fast  dynamics,  polymers  with  greater  mobility 
are  more  susceptible  to  confinement,  especially  below 
the  calorimetric  T . In  contrast,  polymers  with  strongly 
localized  motions  appear  bulk-like,  even  in  highly 
confined  glassy  films.  More  details  of  these  studies  can 
be  found  in  our  upcoming  publications  [3-5]. 
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intimate  Relationship  Between  Structural 
Deformation  and  Properties  of  Carbon  Nanotubes 


Carbon  nanotubes  continue  to  surprise  scientists  with 
their  novel  properties.  Recently  we  have  discovered 
many  intimate  relationships  between  structural  deforma- 
tion and  the  properties  of  single-walled  nanotubes 
(SWNT)  [1,3]  that  could  be  important  in  technological 
applications.  These  observations  were  made  using  state- 
of-the-art  first-principles  total  energy  calculations  within 
the  pseudo-potential  and  the  generalized  gradient  ap- 
proximations (GGA).  Here  we  present  a brief  review  of 
our  findings. 

SWNTs  are  basically  rolled  graphite  sheets  which  are 
characterized  by  two  integers  (tun)  defining  the  rolling 
vector  of  graphite.  Therefore,  the  electronic  properties  of 
SWNTs,  to  first  order,  can  be  deduced  from  that  of 
graphite  by  mapping  the  band  structure  from  a 2D 
hexagonal  lattice  onto  a cylinder.  Such  an  analysis 
indicates  that  (n,n)  armchair  nanotubes  are  always 
metallic  and  (n,0)  zigzag  nanotubes  are  metallic  only  if  n 
is  an  integer  multiple  of  three.  From  our  first-principles 
calculations  we  find  that  these  electronic  properties  can 
be  significantly  modified  by  radial  strain,  which  distorts 
the  circular  cross  section  to  an  elliptical  one.  Figure  1 
shows  that  the  energy  gap  of  an  insulating  SWNT 
decreases  and  eventually  vanishes  at  an  insulator-metal 
transition  with  increasing  applied  elliptical  strain.  The 
density  of  states  at  the  Fermi  level  N(Ef)*  of  a metallized 
SWNT  increases  with  increasing  deformation.  More 
interestingly,  the  elliptical  deformation  necessary  to 
induce  metallicity  is  found  to  be  in  the  elastic  range. 
Therefore,  all  strain  induced  changes  in  electronic  and 
also  in  mechanical  properties  are  reversible. 

Another  significant  effect  of  the  elliptical  distortion  of 
SWNTs  is  the  change  in  the  uniformity  of  their  charge 
distribution.  This,  in  turn,  imposes  changes  in  the 
chemical  reactivity  and  hence  on  the  interaction  of  the 
tube  surface  with  foreign  atoms  and  molecules.  To 
demonstrate  this,  we  calculate  the  binding  energies  of  H 
and  a simple  metal  atom,  Al,  on  SWNTs  as  functions  of 
nanotube  radius  (/?),  type,  and  elliptical  deformation.  We 
find  that  the  binding  energy,  £B,  is  positive  (i.e.  stable) 
and  follows  a remarkably  simple  1 /R  scaling  for  (n,0) 
SWNTs  for  both  H and  Al,  even  though  H prefers  to  sit 


Fig.  1 . Variation  of  the  energy  band  gap  Egap( a)  and  density  of  the 
states  at  the  Fermi  energy  (b)  as  a function  of  elliptical  deforma- 
tion 

on  top  of  C atoms  and  Al  favors  the  hollow  sites  (i.e.  top 
of  hexagons).  Increasing  binding  energy  with  increasing 
curvature  (i.e.  decreasing  R ) suggests  that  EB  can  be 
modified  by  elliptical  deformation. 

Figure  2 shows  the  variation  with  elliptical  distortion 
of  the  binding  energy  E of  a single  hydrogen  atom 
adsorbed  on  the  sharp  and  flat  edges  of  the  (8,0)  surface. 
With  increasing  distortion,  the  binding  energy  of  the 
sharp  site  increases,  while  it  decreases  for  the  fiat  site, 
creating  an  energy  difference  of  ~ 1.1  eV  at  the  distortion 
value  e=  0.3.  This  value  is  44  % of  the  binding  energy  of 
H on  the  undeformed  SWNT. 

The  binding  energy  of  Al  exhibits  a behavior  similar 
to  that  of  H,  despite  H and  Al  favoring  different  sites  on 
the  (8.0)  tube.  Elliptical  deformation  increases  the 
binding  energy  of  the  both  fiat  and  sharp  sites  for  Al, 
given  an  energy  difference  of  0.5  eV  at  the  deformation 
e = 0.3.  It  is  remarkable  that  Al.  which  is  not  bound  to  the 
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Fig.  2.  Binding  energies  of  the  sharp  (top)  and  the  flat  sites  (bot- 
tom) versus  elliptical  deformation  for  a hydrogen  adsorbed  on  a 
(8,0)  SWNT. 

graphite,  can  be  adsorbed  at  the  sharp  site  of  the  distorted 
SWNT  with  a binding  energy  of  1 .8  eV.  Hence,  we 
conclude  that  not  only  band  gap  engineering  but  also 
chemical  reactions  taking  place  on  the  surface  of  a 
SWNT  can  be  engineered  through  radial  deformation. 

Explanation  of  this  remarkable  change  of  the  binding 
energy  with  elliptical  deformation  is  sought  in  the 
electronic  structure  of  SWNTs.  Figure  3 shows  the 
response  of  the  first  conduction  band  to  elliptical  distor- 
tion. It  is  clear  that  the  distortion  disturbs  the  uniformity 
of  the  charge  distribution  of  the  SWNT  and  pushes  the 
chemically  most  active  electrons  from  the  flat  to  the 
sharp  site  of  the  nanotube,  increasing  the  adsorption 
energy  of  the  sharp  site  significantly.  Surprisingly  the 
effect  of  distortion  on  metallic  armchair  (n,n)  nanotubes 
is  found  to  be  very  small  due  to  the  metallic  bands  of  the 
tubes  [2],  This  could  be  important  in  selective 
functionalization  of  nanotubes.  See  Ref.|  1]  for  animation 
of  this  effect. 

We  believe  that  the  tunable  electronic  and  chemical 
properties  reported  here  can  have  important  implications 
for  metal  coverage  and  for  selective  adsorption  and 
desorption  of  atoms  and  molecules  on  carbon  nanotubes. 
This  could  allow  the  fine  tuning  of  the  properties  of 
SWNTs  via  reversible  deformation  and  can  ultimately 
lead  to  a wide  variety  of  technological  applications  such 


Fig.  3.  Contour  plot  of  the  first  conduction  band  of  an  undeformed 
(top)  and  elliptically  deformed  (bottom)  (8,0)  SWNT. 

as  variable  metal-insulator  junctions,  quantum  wells, 
catalysts,  hydrogen  storage  devices,  magnetic  tubes  (by 
absorbing  magnetic  ions),  etc. 
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Serving  the  Scientific  and  Technological 

Communities 


The  role  of  the  NCNR  as  a national  user  facility  has 
expanded  significantly  over  the  past  few  years,  as  the  last 
of  the  instruments  envisioned  in  the  original  cold  neutron 
project  have  become  operational.  The  Disk-Chopper 
Spectrometer,  the  Filter-Analyzer  Neutron  Spectrometer, 
the  High  Flux  Backscattering  Spectrometer,  and  the 
Neutron  Spin-Echo  spectrometer  now  permit  U.S. 
scientists  to  carry  out  neutron  spectroscopy  with  greatly 
improved  resolution  and  sensitivity.  The  thermal  perfect- 
crystal  small  angle  neutron  scattering  instrument 
(USANS)  has  expanded  the  length  scale  available  by  this 
technique  to  104  nm.  User  experiments  show  a steadily 
increasing  diversification  in  subject  area  and  technique, 
enabled  by  the  new  instruments.  We  anticipate  that  the 
trend  will  continue  over  the  next  few  proposal  cycles. 

User  participation  over  the  past  1 5 years  shows 
continuing  growth  (see  Fig.  1 ).  The  NCNR  currently 
accommodates  more  than  half  of  all  neutron  users  in  the 
United  States.  While  the  Spallation  Neutron  Source  is 
being  built  at  Oak  Ridge,  the  NCNR  continues  to  be  the 
nation's  premier  facility  for  providing  neutrons  to  the 
U.S.  research  community. 


In  FY  2001  Research  Participants  were  from 
22  NIST  Divisions  and  Offices 
34  U.S.  Government  Laboratories 
119  U.S.  Universities 
50  U.S.  Industrial  Laboratories 
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Fig.  1.  Research  Participants  at  the  NCNR. 


The  NCNR  User  Program 

Researchers  may  obtain  use  of  NCNR  neutron  beam 
instruments  in  several  ways,  the  most  direct  being 
through  the  formal  proposal  system.  Approximately  twice 
a year,  a Call  for  Proposals  is  issued.  After  a thorough 
review  process  by  external  referees  and  by  the  NCNR 
Program  Advisory  Committee  (PAC),  approved  proposals 
are  allocated  beam  time.  The  PAC  is  a panel  of  distin- 
guished scientists  with  expertise  across  a broad  range  of 
neutron  methods  and  scientific  disciplines.  It  is  the  body 
primarily  responsible  for  proposal  review  and  recom- 
mending user  policies  for  the  NCNR,  working  closely 
with  the  Center’s  Director  and  staff.  Its  current  member- 
ship includes  Sanat  Kumar  (Penn  State  University,  chair), 
Robert  Briber  (University  of  Maryland).  Michael 
Crawford  (DuPont),  Dieter  Schneider  (Brookhaven 
National  Laboratory).  Kenneth  Herwig  (Oak  Ridge 
National  Laboratory),  Yumi  Ijiri  (Oberlin  College), 
Michael  Kent  (Sandia  National  Laboratories,  John 
Tranquada  (Brookhaven  National  Laboratory),  and 
Andrew  Allen  (NIST). 

At  the  recent  PAC  meeting  in  January  2001,  the  PAC 
considered  121  proposals  for  SANS  and  reflectometry,  in 
addition  to  49  for  inelastic  neutron  scattering.  Although 
we  expect  that  both  categories  will  see  increased  user 
demand  in  future  proposal  rounds,  the  latter  area  is  likely 
to  see  more  growth,  since  the  new  inelastic  scattering 
spectrometers  have  added  capabilities  previously  unavail- 
able at  US  neutron  facilities. 

The  Center  for  High  Resolution  Neutron 
Scattering 

Several  NCNR  instruments  are  supported  by  the 
National  Science  Foundation  (NSF)  through  the  Center 
for  High  Resolution  Neutron  Scattering  (CHRNS),  a very 
important  component  of  the  user  program.  Until  the 
present,  CHRNS  instruments  included  a 30  m SANS 
machine,  the  SPINS  triple-axis  spectrometer,  and 
USANS.  Approximately  40  % of  the  instrument  time 
allocated  by  the  PAC  went  to  experiments  carried  out  on 
CHRNS  instruments.  During  the  past  year,  the  NSF 


approved  a proposal  to  provide  support  for  a substantial 
enlargement  of  CHRNS,  enabling  us  to  offer  more  beam 
time  to  users  on  the  three  new  high-resolution  inelastic 
scattering  instruments.  In  the  near  future,  another  SANS 
diffractometer,  the  8 m machine  on  neutron  guide  NG-1, 
which  is  presently  used  primarily  for  NIST  programmatic 
research,  will  be  upgraded  to  a more  powerful  9 m 
instrument  with  a new  detector,  and  made  available  to 
users  through  CHRNS.  The  expanded  suite  of  CHRNS 
instruments  will  provide  enhanced  technical  support  and 
accessibility  to  the  user  community. 

Summer  School  on  Methods  and 
Applications  of  Neutron  Spectroscopy 

The  NCNR's  seventh  annual  Summer  School  on 
neutron  scattering,  held  on  June  18-22,  2001,  focused  on 
high-resolution  neutron  spectroscopy.  Most  of  the  32 
participants  were  graduate  and  postdoctoral  students;  we 
also  had  an  undergraduate  student  and  a professor. 

In  a radical  departure  from  the  format  of  previous 
NCNR  summer  schools,  plenary  lectures  were  presented 
on  the  first  morning  only.  After  lunch  the  students  were 
divided  into  four  groups  of  eight,  and  each  group  per- 
formed an  experiment  using  a different  instrument  (in 
one  case  a pair  of  instruments).  Students  learned  the 
scientific  motivation  for  the  experiment,  how  the  instru- 
ment works,  and  how  to  set  up  measurements.  The 
instruments  collected  data  overnight  and  on  the  following 
morning  the  data  were  analyzed.  After  lunch  on  each  of 
the  second,  third  and  fourth  days  the  groups  rotated  to 
new  experiments.  After  lunch  on  the  final  day,  eight  of 
the  students  gave  short  presentations  describing  an 
experiment. 

Students  used  the  following  instruments:  the  Spin 
Polarized  Inelastic  Neutron  Spectrometer,  the  Disk 
Chopper  Spectrometer,  the  High-Flux  Backscattering 
Spectrometer,  the  Neutron  Spin  Echo  spectrometer,  and 
the  Filter  Analyzer  Neutron  Spectrometer.  Together  these 
instruments  probe  motions  in  materials  over  six  decades 
in  time  (from  approximately  I0‘7  s to  lO13  s)  and  two 
orders  of  magnitude  in  distance  (0.1  nm  to  10  nm).  Test 


Fig.  2.  Seventh  annual  Summer  School  on  Neutron  Scattering 
participants  analyze  results  with  NCNR’s  Seung-Hun  Lee  (bot- 
tom left). 


samples  included  a geometrically  frustrated  antiferroinag- 
net,  a transition  metal  compound  with  organic  ligands,  a 
"symmetric  top"  molecule,  and  a system  of  spherical 
micelles.  The  experiments  illustrated  coherent  scattering, 
incoherent  scattering,  magnetic  scattering,  and  tunneling. 
In  evaluating  the  school  more  than  80  % of  the  partici- 
pants described  it  as  “excellent",  and  generally  high 
ratings  were  received  in  all  areas.  The  lectures  and 
experiment  handouts  have  been  placed  on  the  Web 
( www.ncnr.nist.gov/programs/spectroscopy/SSO  1 / 

SSO l_materials.html.)  A short  report  about  the  summer 
school  will  appear  in  "Neutron  News".  As  in  the  past,  the 
summer  school  was  jointly  sponsored  with  the  National 
Science  Foundation,  which  provided  financial  assistance 
to  many  of  the  university  participants. 

Collaborations 

Direct  collaborations  on  specific  experiments  remain 
a common  way  for  users  to  pursue  their  ideas  using 
NCNR  facilities,  accounting  for  approximately  half  of  the 
number  of  instrument-days.  The  thermal-neutron  triple- 
axis spectrometers  are  mainly  scheduled  in  this  way. 

Most  of  the  time  reserved  for  NIST  on  these  and  all  other 
NCNR  instruments  is  devoted  to  experiments  that  are 
collaborations  with  non-NIST  users.  Collaborative 
research  involving  external  users  and  NIST  scientists 
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often  produces  results  that  could  be  not  obtained  other- 
wise. 

Another  mode  of  access  to  the  NCNR  is  through 
Participating  Research  Teams  (PRTs).  In  this  case,  groups 
of  researchers  from  various  institutions  join  forces  to 
build  and  operate  an  instrument.  Typically,  50  % to  75  % 
of  the  time  on  the  instrument  is  then  reserved  for  the 
PRT.  and  the  remaining  time  is  allocated  to  general  user 
proposals.  For  example,  a PRT  involving  ExxonMobil, 
the  University  of  Minnesota,  and  NIST  cooperates  on  the 
NG-7  30  m SANS  instrument.  Similar  arrangements 
involving  other  PRTs  apply  for  the  horizontal-sample 
reflectometer,  the  high-resolution  powder  diffractometer, 
the  filter-analyzer  spectrometer,  and  the  neutron  spin- 
echo  spectrometer. 

Independent  Programs 

There  are  a number  of  programs  of  long  standing 
located  at  the  NCNR  that  involve  other  parts  of  NIST, 
universities,  industrial  laboratories,  or  other  government 
agencies. 

The  Polymers  Division  of  the  Materials  Science  and 
Engineering  Laboratory  has  two  major  program  elements 
at  the  NCNR.  In  the  first,  the  purpose  is  to  help  the  U.S. 
microelectronics  industry  in  addressing  their  most 
pressing  materials  measurement  and  standards  issues.  In 
today's  integrated  circuits  and  packages  the  feature  size 
on  a chip  is  ever  shrinking,  approaching  250  nm,  while 
the  size  of  a polymer  molecule  is  typically  5 nm  to  10 
nm.  As  feature  size  shrinks,  the  structure  and  properties 
of  interfaces  play  an  increasingly  important  role  in 
controlling  the  properties  of  the  polymer  layers  used  in 
interconnects  and  packages.  NIST  scientists  use  both 
neutron  reflectivity  and  other  neutron  scattering  methods 
to  characterize  polymer/metal  interfaces  with  regard  to 
local  chain  mobility,  moisture  absorption,  glass  transition 
temperature  and  crystalline  structure. 

In  the  second  program  element,  the  objective  is  to 
understand  underlying  principles  of  phase  behavior  and 


phase  separation  kinetics  of  polymer  blends,  both  in  the 
bulk  and  on  surfaces,  in  order  to  help  control  morphology 
and  structure  during  processing.  SANS  and  reflectivity 
measurements  in  equilibrium,  in  transient  conditions,  and 
under  external  fields,  provide  essential  information  for 
general  understanding  as  well  as  for  specific  application 
of  polymer  blend/alloy  systems.  Customers  include 
material  producers  and  users,  ranging  from  chemical, 
rubber,  tire,  and  automotive  companies,  to  small  molding 
and  compounding  companies.  The  focus  of  research  on 
polymeric  materials  includes  commodity,  engineering  and 
specialty  plastic  resins,  elastomers,  coatings,  adhesives, 
films,  foams,  and  fibers. 

The  ExxonMobil  Research  and  Engineering 
Company  is  a member  of  the  Participating  Research 
Team  (PRT)  that  operates,  maintains,  and  conducts 
research  at  the  NG-7  30  m SANS  instrument  and  the 
recently  commissioned  NG-5  Neutron  Spin  Echo  Spec- 
trometer. The  mission  is  to  use  those  instruments,  as  well 
as  other  neutron  scattering  techniques,  in  activities  that 
complement  research  at  ExxonMobil's  main  laboratories 
as  well  as  at  laboratories  of  its  affiliates  around  the 
world.  The  aim  of  these  activities  is  to  deepen  under- 
standing of  the  nature  of  ExxonMobil's  products  and 
processes,  so  as  to  improve  customer  service  and  to 
improve  the  return  on  shareholders'  investment.  Accord- 
ingly, and  taking  full  advantage  of  the  unique  properties 
of  neutrons,  most  of  the  experiments  use  SANS  or  other 
neutron  techniques  to  study  the  structure  and  dynamics  of 
hydrocarbon  materials,  especially  in  the  fields  of  poly- 
mers, complex  fluids,  and  petroleum  mixtures. 
ExxonMobil  regards  its  participation  in  the  NCNR  and 
collaborations  with  NIST  and  other  PRT  members  not 
only  as  an  excellent  investment  for  the  company,  but  also 
as  a good  way  to  contribute  to  the  scientific  health  of  the 
Nation. 

The  Nuclear  Methods  Group  (Analytical  Chemistry 
Division,  Chemical  Science  and  Technology  Laboratory) 
has  as  its  principal  goals  the  development  and  application 
of  nuclear  analytical  techniques  for  the  determination  of 
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elemental  compositions  with  greater  accuracy,  higher 
sensitivity  and  better  selectivity.  A high  level  of  compe- 
tence has  been  developed  in  both  instrumental  and 
radiochemical  neutron  activation  analysis  (INAA  and 
RNAA).  In  addition,  the  group  has  pioneered  the  use  of 
cold  neutron  beams  as  analytical  probes  with  both 
prompt  gamma  activation  analysis  (PGAA)  and  neutron 
depth  profiling  (NDP).  PGAA  measures  the  total  amount 
of  a particular  analyte  present  throughout  a sample  by  the 
analysis  of  the  prompt  gamma-rays  emitted  during 
neutron  capture.  NDP,  on  the  other  hand,  determines 
concentrations  of  several  important  elements  (isotopes) 
as  a function  of  depth  within  the  first  few  micrometers  of 
a surface  by  energy  analysis  of  the  prompt  charged- 
particles  emitted  during  neutron  bombardment.  These 
techniques  (INAA,  RNAA,  PGAA,  and  NDP)  provide  a 
powerful  combination  of  complementary  tools  to  address 
a wide  variety  of  analy  tical  problems  of  great  importance 
in  science  and  technology,  and  are  used  to  help  certify  a 
large  number  of  NIST  Standard  Reference  Materials. 

During  the  past  several  years,  a large  part  of  the 
Group’s  efforts  has  been  directed  towards  the  exploita- 
tion of  the  analytical  applications  of  the  guided  cold- 
neutron  beams  available  at  the  NIST  Center  for  Neutron 
Research.  The  Group's  involvement  has  been  to  design 
and  construct  state-of-the-art  cold  neutron  instruments 
for  both  PGAA  and  NDP  and  provide  facilities  and 
measurements  for  outside  users,  while  retaining  and 
utilizing  our  existing  expertise  in  INAA  and  RNAA. 

The  Center  for  Food  Safety  and  Applied  Nutrition. 
U.S.  Food  and  Drug  Administration  (FDA),  directs  and 
maintains  a neutron  activation  analysis  (NAA)  facility  at 
the  NCNR.  This  facility  provides  agency-wide  analytical 
support  for  special  investigations  and  applications 
research,  complementing  other  analytical  techniques  used 
at  FDA  with  instrumental,  neutron-capture  prompt- 
gamma,  and  radiochemical  NAA  procedures,  radioiso- 
tope x-ray  fluorescence  spectrometry  (RXRFS),  and 
low-level  gamma-ray  detection.  This  combination  of 
analytical  techniques  enables  diverse  multi-element  and 


radiological  information  to  be  obtained  for  foods  and 
related  materials.  The  NAA  facility  supports  agency 
quality  assurance  programs  by  developing  in-house 
reference  materials,  by  characterizing  food-related 
reference  materials  with  NIST  and  other  agencies,  and  by 
verifying  analyses  for  FDA's  Total  Diet  Study  Program. 
Other  studies  include  the  development  of  RXRFS 
methods  for  screening  foodware  for  the  presence  of  Pb. 
Cd  and  other  potentially  toxic  elements,  use  of  instru- 
mental NAA  to  investigate  bromate  residues  in  bread 
products,  and  use  of  prompt-gamma  NAA  to  investigate 
boron  nutrition  and  its  relation  to  bone  strength. 

The  Neutron  Interactions  and  Dosimetry  Group 
(Physics  Laboratory)  provides  measurement  services, 
standards,  and  fundamental  research  in  support  of  NIST's 
mission  as  it  relates  to  neutron  technology  and  neutron 
physics.  The  national  and  industrial  interests  served 
include  scientific  instrument  calibration,  electric  power 
production,  radiation  protection,  defense  nuclear  energy 
systems,  radiation  therapy,  neutron  radiography,  and 
magnetic  resonance  imaging. 

The  Group's  activities  may  be  represented  as  three 
major  activities.  The  first  is  Fundamental  Neutron 
Physics  including  magnetic  trapping  of  ultracold  neu- 
trons, operation  of  a neutron  interferometry  and  optics 
facility,  development  of  neutron  spin  filters  based  on 
laser  polarization  of 'He,  measurement  of  the  beta  decay 
lifetime  of  the  neutron,  and  investigations  of  other 
coupling  constants  and  symmetries  of  the  weak  interac- 
tion. This  project  involves  a large  number  of  collabora- 
tors from  universities  and  national  laboratories. 

The  second  is  Standard  Neutron  Fields  and  Applica- 
tions utilizing  both  thermal  and  fast  neutron  fields  for 
materials  dosimetry  in  nuclear  reactor  applications  and 
for  personnel  dosimetry  in  radiation  protection.  These 
neutron  fields  include  thermal  neutron  beams,  “white" 
and  monochromatic  cold  neutron  beams,  a thermal- 
neutron-induced  235U  fission  neutron  field,  and  252Cf 
fission  neutron  fields,  both  moderated  and  unmoderated. 


The  third  is  Neutron  Cross  Section  Standards  includ- 
ing experimental  advancement  of  the  accuracy  of  neutron 
cross  section  standards,  as  well  as  evaluation,  compila- 
tion and  dissemination  of  these  standards. 

Several  universities  have  also  established  long  term 
programs  at  the  NCNR.  The  University  of  Maryland  is 
heavily  involved  in  the  use  of  the  NCNR.  and  maintains 
several  researchers  at  the  facility.  Johns  Hopkins 
University  participates  in  research  programs  in  solid- 
state  physics  and  in  instrument  development  at  the 
NCNR.  The  University  of  Pennsylvania  is  working  to 
help  develop  biological  applications  of  neutron  scatter- 
ing. It  is  also  participating  in  the  second  stage  construc- 
tion of  the  filter  analyzer  neutron  spectrometer,  along 
with  the  University  of  California  at  Santa  Barbara, 
DuPont,  Hughes,  and  Allied  Signal.  The  University  of 
Minnesota  participates  in  two  PRTs,  the  NG-7  30  m 
SANS  and  the  NG-7  reflectometer.  The  University  of 
Massachusetts  also  participates  in  the  latter  PRT. 


Reactor  Operation  and  Engineering 


The  reactor  operated  for  240  full  power  (20  MW) 
days  during  the  past  year,  or  1 00  % of  the  scheduled 
time.  Routinely,  the  reactor  is  scheduled  to  operate  a 
seven  week  cycle,  seven  times  a year.  Each  operating 
cycle  includes  38  days  on-line  and  1 1 days  shutdown  for 
refueling,  routine  maintenance,  and  surveillance  tests. 
(Figure  I shows  the  reactor  operations  staff  adjusting  the 
reactor  power  regulatory  rod).  This  year  four  days  in 
March  were  required  for  maintenance  on  the  Thermal 
Column  Cooling  System.  The  reactor  was  shutdown  on 
August  26,  2001  for  a planned  maintenance  period  of 
three  months.  During  this  time  the  cold  source  will  be 
replaced,  building  transformers  and  switchgear  will  be 
replaced  to  provide  additional  electrical  capacity,  and  the 
new  cooling  tower  will  undergo  final  hookup  and  testing. 
(Figure  2 shows  the  new  cooling  tower  in  the  final  stages 
of  construction.  For  comparison  the  old  cooling  tower  is 
seen  on  the  left). 

Analysis  has  been  completed  to  verify  that  the 
thermal-hydraulic  characteristics  of  the  NBSR,  as  was 
predicted  using  correlations  and  simple  models  available 
when  the  licensed  power  level  was  raised  to  20  MW,  are 
conservative  relative  to  results  obtained  using  newer 
correlations  that  better  represent  research  reactor  charac- 
teristics. Calculations  show  that  the  heat  transfer  coeffi- 
cient, onset  of  How  instability,  and  critical  heat  flux  are 
all  conservative.  In  addition,  the  NBSR  also  has  an 
adequate  margin  before  the  onset  of  nuclear  boiling. 

The  above  results  are  for  extreme  yet  still  conceiv- 
able conditions  of  operation.  The  Energy  Sciences  and 
Technology  Department  at  the  Brookhaven  National 
Laboratory  has  been  tasked  to  perform  transient  analysis 
with  state-of-the  art  methods.  Monte  Carlo  neutronics 
calculations  will  be  used  to  accurately  predict  parameters 
that  are  used  in  the  transient  analysis. 

Work  is  also  continuing  in  the  preparation  of  an 
updated  license  renewal  submittal  to  the  NRC.  In  con- 
junction with  this  effort,  an  in-service  inspection  of 
reactor  internals  and  ultrasonic  testing  of  the  primary 
cooling  system  are  planned  for  the  summer/fall  shut- 
down. 


Fig.  1.  Bill  Mueller  and  Greg  Heller  of  the  reactor  operations  staff 
adjust  the  reactor  power  regulatory  rod. 


Fig.  2.  The  new  plume-abated  cooling  tower  (right)  nears 
completion.  It  will  replace  the  lower  capacity  and  less  efficient 
system  that  had  no  abatement  ability  (left). 


Instrumentation  Developments 


The  Advanced  Cold  Neutron  Source 

After  6 years  of  safe  and  reliable  operation,  the  liquid 
hydrogen  cold  source  cryostat  at  the  NIST  reactor  is 
being  replaced  with  a new  advanced  design.  Based  on 
extensive  operating  experience  and  sophisticated  Monte 
Carlo  calculations  of  the  cold  source  and  NIST  reactor 
core,  the  new  design  improves  the  neutron  coupling 
between  the  reactor  fuel  and  moderator  and  the  cold 
source  cryostat.  The  advanced  source,  (called  Unit  2), 
will  differ  from  the  original  in  many  key  respects,  as  can 
be  seen  in  Fig.  I . The  most  important  change  is  the 
addition  of  60  liters  of  D.O.  partially  surrounding  the 
moderator  chamber,  reducing  the  moderator  void  volume 
and  increasing  the  thermal  neutron  flux  in  the  cryostat 
region  by  about  40  %. 


Fig.  1.  Comparison  between  the  original  (top)  and  new  designs 
for  the  liquid  hydrogen  cold  source  cryostat  at  NIST.  Heavy  wa- 
ter volume  is  indicated  in  blue,  liquid  hydrogen  in  orange,  and 

hydrogen  vapor  in  yellow. 


The  new  cryostat  is  shaped  like  an  ellipsoidal  shell, 
rather  than  a spherical  one.  This  results  in  a smaller 
volume  so  that  more  D.O  can  be  introduced  in  the 
cryostat  assembly.  The  center  of  the  inner  ellipsoid  is 
offset,  so  the  liquid  hydrogen  shell  is  30  mm  thick  near 
the  reactor  core.  This  change  results  in  an  additional 
50  % increase  in  flux.  Finally,  in  contrast  to  the  original 
design,  the  inner  ellipsoid  will  be  evacuated  through  a 
small  port,  removing  most  of  the  hydrogen  vapor.  This 
change  will  add  another  1 5 % to  20  % gain.  Based  on 
these  changes,  users  can  expect  an  overall  intensity 


Fig.  2.  False  color  representation  of  stresses  from  finite  element 
calculations  of  the  helium  containment  jacket  with  an  internal 
pressure  of  82  atmospheres. 

increase  of  50  % to  80  % between  2 A and  20  A with  this 
design. 

The  complex  geometry  of  the  advanced  source  posed 
many  engineering  challenges.  The  moderator  chamber  is 
surrounded  by  a vacuum  vessel,  which  is  surrounded  by  a 
helium  containment  vessel,  strong  enough  to  withstand 
the  design  basis  accidental  detonation  of  liquid  hy  drogen 
and  solid  oxygen.  These  vessels  determine  the  extent  of 
the  D.O  volume,  so  they  were  designed  to  partially 
surround  the  moderator  chamber,  allowing  an  opening 
just  large  enough  to  fully  illuminate  the  neutron  guides. 


This  “hour-glass”  shape  created  regions  of  very  high 
mechanical  stress,  analyzed  and  ultimately  mitigated  with 
the  aid  of  finite  element  analysis  (see  Fig.  2).  During  the 
three  month  period  beginning  at  the  end  of  August  200 1 . 
we  are  removing  the  old  cryostat  and  installing  Unit  2. 

Development  of  an  Air-pad  and  Floor 
System  for  New  Neutron  Spectrometers 

The  NCNR  is  currently  developing  several  new 
neutron  triple  axis  spectrometers  with  a key  design 
feature  that  allows  interchangeability  of  the  back-end 
detector/analyzer  systems.  Exchanging  the  back-end 
system  on  these  instruments  requires  a universal  connec- 
tion and  positioning  system  between  it  and  the  sample 
position  of  the  instrument.  This  interface  must  be  capable 
of  achieving  an  angular  accuracy  of  0. 1 degrees,  an 
angular  precision  of  0.01  degrees  and  a vertical  posi- 
tional accuracy  of  less  than  a millimeter. 

The  design  approach  adopted  uses  a commercial  air- 
pad  system  with  a precision-leveled  floor.  After  consider- 
ing several  different  types  of  floors  we  chose  to  use 
epoxy  as  our  sub-floor  with  precision  cut  and  anodized 
aluminum  plates  as  our  contact  surface.  Epoxy  will  self- 
level, fill  in  holes,  and  mask  over  any  imperfections  in 
the  concrete  floor  and  results  in  a highly  smooth  and 
level  surface  to  build  on.  The  epoxy  is  also  low  cost. 


Fig.  3.  Michael  Murbach  and  Colin  Wrenn  conduct  tests  on  the 
prototype  floor  system  for  a triple  axis  spectrometer. 


relatively  easy  to  use,  and  provides  great  flexibility  for 
repair  or  maintenance.  The  aluminum  plates  provide 
protection  to  the  epoxy  from  cryogenic  spills,  increased 
wear  resistance,  and  an  easily  removable/reusable  floor. 

We  have  developed  and  tested  a prototype  single  air- 
pad  system  that  has  successfully  demonstrated  the 
construction  techniques  for  the  floor  system  and  the 
ability  of  the  air  pad  positioning  system  to  meet  the 
accuracy  and  repeatability  requirements  for  a 1 ton  test 
load.  Currently,  a three-pad  prototype  system  is  under 
development  which  will  allow  us  to  test  control  system 
concepts  needed  in  a final  design. 

Development  of  an  Intercalated  Graphite 
Monochromator 

As  part  of  an  ongoing  collaboration  between  NIST 
and  Harvard  University  to  measure  the  neutron  lifetime 
using  magnetically  trapped  ultracold  neutrons  (UCNs),  a 
new  monochromator  has  been  developed.  In  the  lifetime 
experiment,  ultracold  neutrons  (UCN:  neutrons  with 
energies  ~ 100  neV)  are  produced  in  a supertluid  4He 
bath  inside  the  trapping  apparatus  using  the 
“superthermal  process.”  A beam  of  cold  neutrons  is 
introduced  into  the  helium,  in  which  neutrons  with 
wavelength  near  0.89  nm  scatter  to  UCN  energies  by 
creation  of  a single  phonon.  Only  neutrons  within  a 
narrow  wavelength  band  (±  0.01  nm  for  our  current 
magnetic  trap)  participate  in  this  single-phonon  produc- 
tion mechanism.  Since  the  wavelength  band  for  UCN 
production  occupies  only  a small  fraction  of  the  neutron 
spectrum  produced  by  the  cold  source  (see  Fig.  4),  a 
significant  signal-to-background  improvement  is  possible 
with  appropriate  wavelength  filtering  of  neutrons  that 
enter  the  apparatus.  Bragg  reflection  from  a crystalline 
monochromator  with  a lattice  spacing,  t/,  can  be  used  to 
select  a narrow  wavelength  band  for  neutrons  with 
wavelengths  up  to  2d.  Unfortunately,  widely  used 
materials  for  neutron  monochromatization,  such  as 
pyrolytic  graphite,  have  insufficiently  large  J-spacings 
for  reflecting  0.89  nm  neutrons.  One  attractive  class  of 
materials  for  producing  long  wavelength  monochroma- 
tors is  the  intercalated  graphites.  Intercalation  is  the 
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process  of  inserting  foreign  atoms  or  molecules  between 
the  layers  of  the  graphite,  and  results  in  a crystal  with  a 
larger  d spacing  Intercalated  graphites  are  formed  in 
“stages”,  where  the  stage  number  defines  the  number  of 
planes  of  graphite  between  layers  of  intercalant  atoms. 

We  have  designed,  fabricated,  and  tested  a new 
0.89  nin  monochromator  using  stage-2  potassium- 
intercalated  graphite.  Nine  pieces  of  intercalated  graphite 


Fig.  4.  Neutron  wavelength  spectra  were  measured  using  a time- 
of-flight  technique.  Only  1 % of  the  neutrons  emerging  from  guide 
NG-6  (red)  lie  within  the  wavelength  band  (between  the  orange 
vertical  bars)  available  for  production  of  UCN.  The  spectra  of 
neutrons  transmitted  (blue)  and  reflected  (green)  by  the  mono- 
chromator are  also  shown. 

are  tiled  to  form  the  monochromator,  which  has  a total 
size  of  6 cm  x 15  cm.  Each  tile  is  produced  by  heating 
potassium  and  graphite  to  two  different  temperatures  in 
an  evacuated  glass  cell.  Performance  of  the  assembled 
monochromator  was  determined  using  time-of-flight 
techniques  (see  Fig.  4).  The  monochromator  reflects 
more  than  80  % of  the  neutrons  incident  upon  it  in  the 
wavelength  band  relevant  for  UCN  production.  However, 
only  3 % of  all  incident  neutrons  are  reflected.  When 
combined  with  additional  techniques  that  filter  out 
neutrons  from  higher  order  Bragg  reflections,  the  mono- 
chromator results  in  an  improvement  of  the  ratio  of 
useful  to  total  neutrons  of  almost  two  orders  of  magni- 
tude while  sacrificing  few  trapped  neutrons.  During  the 
summer  of  2001,  this  new  monochromator  was  used  in 
conjunction  with  the  trapping  apparatus  in  a temporary 


setup.  Presently,  we  are  permanently  installing  the 
monochromator  at  the  fundamental  physics  station  at 
NG-6.  This  installation  will  place  the  monochromator 
much  closer  to  the  end  of  the  neutron  guide  and  will 
result  in  an  increase  by  an  additional  factor  of  2 in  the 
number  of  trapped  UCNs. 

NCNR  Computer  and  Software 
Environment 

The  NCNR's  computerized  data  acquisition  and 
reduction  environment  encompasses  all  stages  of  re- 
search: experiment  preparation,  data  collection  and 
analysis  of  results.  NCNR  provides  site-specific  tools  that 
augment  the  tools  in  researchers’  own  scientific  arsenal. 
Some  specific  areas  of  recent  activity  are  described 
below: 

• NCNR  is  providing  simulation  tools  that  model 
the  behavior  of  our  instruments,  to  evaluate  the  feasibility 
of  proposed  experiments.  They  can  be  used  by  prospec- 
tive users,  especially  those  who  aren’t  already  familiar 
with  neutron  scattering  experiments,  to  explore  'what-if' 
scenarios  with  different  sample  and  instrument  configura- 
tions. Currently,  we  have  simulations  of  small  angle 
scattering  experiments,  with  reflectivity  and  stress 
measurement  simulations  being  developed. 

• Researchers  who  plan  to  use  NCNR  instrumenta- 
tion must  submit  their  proposals  for  review.  We  have 
created  a web-based  proposal  submission  system  that 
tracks  user  proposals  from  submission  to  review  and 
ranking,  to  the  final  instrument  time  allocation. 

• All  instruments  available  at  NCNR  are  com- 
puter-controlled. The  control  software  running  on  these 
computers  must  satisfy  two  seemingly  contradictory 
needs:  simplicity,  for  straightforward  use,  and  sophistica- 
tion, for  power  uses.  New,  more  complex  instruments, 
such  as  the  spin  echo  spectrometer,  require  carefully 
designed  graphical  user  interfaces  (GUIs)  that  hide 
unnecessary  detail,  while  still  providing  access  to  all 
functions,  if  needed.  Since  the  NCNR  instruments  run 
around  the  clock,  the  data  acquisition  functions  must  be 
capable  of  automation  and  unattended  operation.  To  this 
end,  the  experiment  control  software  is  scriptable.  in 
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Fig.  5.  A screenshot  of  several  NCNR  applications.  Clockwise  ment,  peak  analysis  application,  Time-of-flight  detector  selec- 
from  top  left  corner:  the  web-based  simulator  of  a SANS  experi-  tion,  and  a Neutron  Spin  Echo  instrument  control  panel. 


addition  to  the  GUI  interface.  The  scriptable,  command- 
based  interface  fulfills  “power-user”  needs,  while  a 
simplified  GUI  interface  tends  to  be  easier  to  learn  and 
use. 

• Rapid  post-acquisition  data  reduction  is  impor- 
tant to  guide  the  ongoing  experiments,  and.  later,  to  aid  in 
the  processing  of  data  to  obtain  the  scientific  results. 
Many  data  reduction  procedures  are  common  to  all  the 
instruments.  The  Data  Visualization  and  Analysis  Envi- 
ronment (DAVE)  is  a set  of  tools  commonly  used  by 
NCNR  scientific  community;  these  tools  include  pro- 


grams to  manipulate  and  compare  data,  extract  param- 
eters of  theoretical  models  by  data  fitting,  etc. 

• NCNR  maintains  a sophisticated  computer 
infrastructure  including  a modern  network  installation, 
and  an  ensemble  of  nearly  400  computers  and  other  IT 
equipment,  for  scientific  and  office  computing,  ranging 
from  email  servers  to  supercomputer-style  parallel 
computing  clusters. 

• NCNR  participates  in  inter-laboratory  collabora- 
tions developing  software  for  specialized  domains. 


Perhaps  the  best  known  of  these  is  the  GSAS  software 
for  crystallographic  data  reduction.  Another  ongoing 
NCNR  collaboration  is  in  international  scattering  data 
exchange  format,  NeXuS. 

Sample  Environment  Team 

The  sample  environment  team  expanded  by  two  new 
members  this  year,  bringing  the  total  number  to  four 
employees  that  are  actively  working  to  ensure  that  the 
needs  of  researchers  working  at  the  NCNR  are  met.  This 
team  works  to  improve  the  infrastructure  sample  environ- 
ment equipment  and  to  support  neutrons  scattering 
instruments  at  the  NCNR.  A principal  objective  is  to 
provide  safe  and  easy  to  use  equipment  to  the  research 
staff.  Equipment  development  activities  are  done  in 
collaboration  with  the  engineering  and  technical  staff, 
with  input  from  both  staff  and  guest  scientists.  As  an 
example,  a new  superconducting  magnet/dilution  refrig- 
erator system  needed  modifications  for  use  on  the  Disk 
Chopper  Spectrometer  (DCS)  for  an  experiment.  The 
impracticality  of  separating  the  cryostat  from  its  support 
pumps  and  electronics  led  to  the  decision  to  lift  the  entire 
system  onto  the  instrument  for  use.  A new  adjustable 
lifting  frame  was  design  and  manufactured  to  lift  the 
various  components  of  the  magnet  (see  Figure  6). 


Another  experiment  required  the  development  of  a 
system  capable  of  going  to  temperatures  as  high  as  675  K 
with  a high  voltage  applied  to  the  sample.  The  design 
goal  was  to  be  able  to  achieve  a 5 kV  potential  difference 
across  the  sample.  After  continued  development  work, 
we  were  able  to  apply  a potential  difference  of  8 kV  to 
the  sample  while  at  the  maximum  temperature  of  675  K. 
In  addition  to  the  human  safety  features  designed  in  from 
the  beginning,  a circuit  was  developed  which  prevents 
the  continued  arcing  that  can  occur  if  the  sample  breaks 
down  during  the  course  of  the  experiment. 


Fig.  6.  Superconducting  magnet/dilution  refrigerator  as  it  is  be- 
ing lifted  onto  the  DCS  spectrometer.  The  adjustable  lifting  frame 
compensates  for  the  lifting  point  shifting  in  case  liquid  cryo- 
gens  are  added. 
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Instruments  and  Contacts 


High  resolution  powder  diffractometer  (BT-1) 

J.  K.  Stalick,  (301)  975-6223,  judy.stalick@nist.gov 

B.  H.  Toby,  (301 ) 975-4297,  brian.toby@nist.gov 

DARTS,  Residual  stress  and  texture  diffractometer  (BT-8) 

H.  J.  Prask,  (301 ) 975-6226,  hank@nist.gov 

T.  Gnaeupel-Herold,  (301 ) 975-5380,  thomas.gnaeupel- 

herold@nist.gov 

30-m  SANS  instrument  (NG-7) 

J.  C.  Schulz,  (301 ) 975-6469,  jamie.schulz@nist.gov 

C.  J.  Glinka,  (301)  975-6242,  cglinka@nist.gov 

J.  G.  Barker,  (301)  975-6732,  john.barker@nist.gov 

30-m  SANS  instrument  (NG-3)  (NIST/NSF-CHRNS) 

B.  Hanunouda,  (301 ) 975-3961,  hammouda@nist.gov 
S.  R.  Kline,  (301)  975-6243,  steven.kline@nist.gov 

D.  Ho,  (301 ) 975-6422,  derek.ho@nist.gov 

8-ni  SANS  instrument  (NG-I ) 

D.  Ho,  (301 ) 975-6422,  derek.ho@nist.gov 

C.  J.  Glinka,  (301)  975-6242,  cglinka@nist.gov 

J.  G.  Barker,  (301)  975-6732,  john.barker@nist.gov 

USANS,  Perfect  crystal  SANS  (BT-5)  (NIST/NSF- 
CHRNS) 

J.  G.  Barker,  (301 ) 975-6732,  john.barker@nist.gov 
C.  J.  Glinka,  (301 ) 975-6242,  cglinka@nist.gov 

Cold  neutron  reflectometer-vertical  sample-polarized 
beam  option  (NC-1) 

C.  F.  Majkrzak,  (30 1 ) 975-525 1 , cmajkrzak@nist.gov 

J.  A.  Dura,  (301)  975-625 1 , jdura@nist.gov 

Cold  neutron  reflectometer-horizontal  sample  (NG-7) 

S.  K.  Satija,  (301)  975-5250,  satija@nist.gov 

K.  W.  Shin,  (301)  975-8829,  kwshin@nist.gov 

Triple-axis  polarized-beam  spectrometer  (BT-2) 

J.  W.  Lynn,  (301 ) 975-6246,  jeff.lynn@nist.gov 

R.W.  Erwin,  (301)975-6245,  rerwin@nist.gov 

Triple-axis  fixed  incident  energy  spectrometer  (BT-7) 

J.  W.  Lynn,  (301 ) 975-6246,  jeff.lynn@nist.gov 
R.W.  Erwin,  (301  )975-6245,  rerwin@nist.gov 

Triple-axis  spectrometer  (BT-9) 

R.  W.  Erwin,  (301)  975-6245,  rerwin@nist.gov 

P.  M.  Gehring,  (301 ) 975-3946,  pgehring@nist.gov 

SPINS,  Spin-polarized  triple-axis  spectrometer  (NG-5) 
(NIST/NSF-CHRNS) 

S.  Park,  (301 ) 975-8369,  sungil.park@nist.gov 

P.  M.  Gehring,  (301 ) 975-3946,  pgehring@nist.gov 
S.-H.  Lee,  (301 ) 975-4257,  seung-hun.lee@nist.gov 
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ill 


FANS,  Filter-analyzer  neutron  spectrometer  (BT-4) 

T.  J.  Udovic,  (301)  975-6241,  udovic@nist.gov 

C.  M.  Brown,  (301)  975-5134,  craig.brown@nist.gov 

D.  A.  Neumann,  (301 ) 975-5252,  dan@nist.gov 


FCS,  Fermi-chopper  time-of-flight  spectrometer  (NG-6) 

C.  M.  Brown,  (301 ) 975-5134,  craig.brown@nist.gov 
T.  J.  Udovic,  (301)  975-6241,  udovic@nist.gov 


DCS,  Disk-chopper  time-of-flight  spectrometer  (NG-4) 
(NIST/NSF-CHRNS) 

J.  R.  D.  Copley,  (301 ) 975-5133,  jcopley@nist.gov 
J.  C.  Cook,  (301 ) 975-6403,  jeremy.cook@nist.gov 

HFBS,  High-llux  backscattering  spectrometer  (NG-2) 
(NIST/NSF-CHRNS) 

R.  M.  Dimeo,  (301)  975-8135,  robert.dimeo@nist.gov 
Z.  Chowdhuri,  (301 ) 975-4404,  zema.chowdhuri@nist.gov 
D.  A.  Neumann,  (301)  975-5252,  dan@nist.gov 

NSE,  Neutron  spin  echo  spectrometer  (NG-5)  (NIST/ 
NSF-CHRNS) 

N.  S.  Rosov,  (301 ) 975-5254,  nrosov@nist.gov 

Prompt-gamma  neutron  activation  analysis  (NG-7) 

R.  M.  Lindstrom,  (30 1 ) 975-628 1 , 
dick.lindstrom@nist.gov 
R.  L.  Paul,  (301 ) 975-6287,  rpaul@nist.gov 


Other  activation  analysis  facilities 

R.  R.  Greenberg,  (301 ) 975-6285,  rgreenberg@nist.gov 

Cold  neutron  depth  profiling  (NG-0) 

G.  Lamaze,  (301 ) 975-6202,  lamaze@nist.gov 

Instrument  development  station  (NG-0) 

D.  F.  R.  Mildner,  (301 ) 975-6366,  mildner@nist.gov 

H.  H.  Chen,  (301 ) 975-3782,  chenmayer@nist.gov 

Neutron  interferometer  (NG-7) 

M.  Arif,  (301 ) 975-6303,  muhammad.arif@nist.gov 

Fundamental  neutron  physics  station  (NG-6) 

M.  S.  Dewey,  (301 ) 975-4843,  mdewey@nist.gov 

Theory  and  modeling 

N.  F.  Berk,  (301 ) 975-6224,  nfb@nist.gov 

T.  Yildirim,  (301 ) 975-6228,  taner@nist.gov 

Sample  environment 

D.  C.  Dender,  (301)  975-6225.  dender@nist.gov 
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